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conditions  has  been  produced  detailing  the  conditions  under  which  useful 
vibrational  population  inversions  may  be  expected  in  the  molecules  produced 
from  such  heterogeneous  catalysis. 
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INTRODUCTION 


Contract  F49620-77-C-0004  was  begun  on  October  1,  1976  to  carry 
out  a "Computational  Study  of  Chemical  Reaction  Dynamics  at  the  Gas-Solid 
Interface".  The  chief  goal  of  the  project  was  to  ascertain  the  conditions 
under  which  a heterogeneously  catalyzed  chemical  reaction  would  channel 
useful  amounts  of  reaction  energy  into  vibrational  modes  of  the  products. 
Possibilities  for  novel  chemical  lasers  were  forseen  should  such  heterogeneous 
up-conversion  be  thoroughly  understood.  A three-year  program  was  envisioned 
to  elucidate  by  means  of  detailed  calculations  what  features  would  be 
necessary  for  efficient  up-conversions.  This  goal  was  realized  in  two  years 
due  to  unexpected  progress  in  developing  computational  procedures  to  model 
surface  dynamics  and  the  possibilities  of  gas-surface  energy  transfer. 
Therefore,  this  report  marks  a natural  transition  point  in  the  program.  In 
order  to  have  a more  coherent  presentation,  we  will  summarize  our  overall 
results  and  conclusions  since  the  inception  of  the  project,  indicating  sepa- 
rate second  year  accomplishments  as  part  of  the  unified  discussion  of  two- 
year  program  accomplishments. 
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IL PROJECT  AIMS 

When  chemical  reactions  occur  in  the  gas  phase  with  release  of 

energy,  this  energy  may  appear  preferentially  in  vibration  of  the  product 
(1 1* 

molecule.  ' These  vibrational ly  excited  product  molecules  may  be  used 

to  produce  lasing  action  in  the  infrared  region  of  the  spectrum. Vibra- 

tionally  excited  HF  or  DF  can  be  produced  by  a variety  of  gas  phase  chemical 
(2) 

reactions.'  ' The  large  amount  of  reaction  energy,  preference  for  vibrational 
excitation  of  the  product  molecules,  and  relatively  inefficient  quenching 
with  typical  diluent  gases  have  produced  large  energy  storage  and  high  powers 
in  DF/HF  chemical  lasers.  Favorable  transmission  characteristics  through 
the  atmosphere,  and  relatively  compact  size,  are  two  features  that  have  made 
the  DF  chemical  laser  probably  the  leading  candidate  in  the  DoD  high  energy 
laser  weapons  program  (although  serious  challenges  from  pulsed  CO2  lasers 
and  others  are  emerging) A DF  chemical  laser  has  reportedly  destroyed  a 
high  speed  antitank  missile  in  flight,  encouraging  proponents  of  high  energy 
lasers  as  weapons  systems.'  ' In  summary,  chemical  lasers  seem  to  be  an 
important  area  of  DoD  concern  expected  to  lead  to  prototype  weapons  systems 
in  the  early  to  mid  1980's.  Research  on  molecular  effects  underlying  chemi- 
cally pumped  lasers  and  novel  ways  to  achieve  such  chemically  pumped  vibra- 
tional population  inversions  supports  the  long-term  development  of  this  area. 
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One  chief  goal  in  this  project  is  to  investigate  via  computations 
a novel  way  to  achieve  population  inversions  in  chemically  reactive  systems. 

The  early  work  of  Polanyi^®^  on  vibrational  population  inversions  in  exo- 
thermic chemical  reactions  in  gases  led  to  the  development  of  the  first  chemical 
laser  by  Kasper  and  Pimental.^^^  These  systems  all  involve  gaseous  reactants 
and  products.  The  intervening  18  years  has  witnessed  a large  growth  in  the 
field  of  gas  phase  chemical  dynamics.  Today,  when  gas  phase  molecules  react 

and  liberate  energy,  it  is  reasonably  well  understood  the  conditions  under  i i 

which  considerable  vibrational  population  inversion  in  the  products  will  be  i 

(8) 

obtained.'  ' As  a result,  interest  in  obtaining  "propensity  rules",  (e.g.,  ' ; 

which  reactions  have  propensities  for  converting  reaction  energy  into  vibra- 
tional  energy)  has  waned  as  far  as  gas  phase  bimolecular  reactions  are  con- 
cerned. I 

However,  much  practical  chemistry  is  performed  with  a heterogeneous  j 

I 

catalyst.  Many  commercially  desirable  reactions  are  inefficient,  too  slow, 

or  produce  undesirable  products.  Carefully  chosen  solids  are  frequently  used  ( 

to  promote  such  reactions  in  economically  desirable  ways.  The  automobile  , 

catalytic  converter  is  one  example  in  which  gas  phase  chemical  reactions  are 

t - 

adjusted  to  yield  desirable  products.  The  present  research  program  represents  / ^ 

I 

the  second  year  of  a research  program  designed  to  elucidate  whether,  and  under  j 

what  conditions,  such  heterogeneous  catalyzed  chemical  reactions  could  lead  I ' 

I K j 

^ J 

to  useful  population  inversion.  Therefore,  by  a combination  of  theoretical 

techniques  from  materials  science,  chemical  dynamics  involving  many-particle  | 

systems,  semi -empirical  quantum  mechanics  and  non-equilibrium  statistical 

mechanics  (Appendix  C),  we  model  the  dynamics  of  reactions  occurring  on  solid  ' 

surfaces.  By  analogy  with  gas  phase  chemistry,  we  expect  to  elucidate  the 
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conditions  under  which  useful  vibrational  population  inversion  may  be  obtained. 
As  discussed  below,  we  believe  significant  progress  has  been  made  in  deter- 
mining necessary  conditions  for  population  inversion  to  occur.  The  next  step 
towards  achieving  a novel  laser  pumping  mechanism  is  for  careful  experiments 
to  be  carried  out  testing  and  refining  the  hypothesis  proposed  here.  One  such 
experimental  program  is  underway  at  Yale.'  ' The  measurement  of  vibrational 
state  distributions  for  comparison  with  our  current  hypothesis  is  one  way  this 
experimental  program  is  likely  to  evolve,  having  benefited  considerably  from 
our  computational  results  as  reported  here  and  in  the  open  literature. 
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III.  PROJECT  ACCOMPLISHMENTS 

A.  Sutimarv 


We  have  performed  a series  of  computer  simulations  that  indicate 
vibrational  excitation  should  occur  for  certain  classes  of  heterogeneous 
reactions.  The  necessary  restrictions  are: 

1.  the  reactants  should  dissociatively  chemisorb  on  the 

surface 

2.  the  product  molecules  should  not  dissociately  chemisorb 

3.  the  masses  of  atoms  going  into  the  product  molecule 

should  be  comparable 

4.  the  energy  that  is  released  should  be  released  early  in 

the  recombination  step 

5.  the  solid  should  be  such  that  energy  transfer  between  the 

gas  and  solid  is  minimized. 

For  these  conditions,  we  have  estimated  that  "67%  of  the  available  reaction 

energy  could  be  channeled  into  vibration  and,  of  this  vibrational  energy, 

about  60%  appears  in  vibrational  levels  V > 2.  The  total  amount  of  energy 

appearing  as  vibration  compares  favorably  with  the  HF  laser  (also  about  67%)^^^ 

but  the  HF  laser  can  be  made  to  produce  more  of  the  energy  in  V ^ 2 (approxi- 
f 3l 

mately  85%).'  ' While  our  model  calculations  have  not  produced  vibrational  popu- 
lations clearly  superior  to  the  existing  (and  highly  optimized)  HF  laser,  use- 
ful inversions  do  seem  to  be  possible. 


2556 
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B.  Details 


A complete  model  of  heterogeneous  reaction  dynamics  involves 
several  steps: 

1.  The  reactants  collide  with  the  solid,  dissociatively 

chemisorb  and  transfer  a portion  of  their  kinetic 
energy  and  energy  of  chemisorption  to  the  solid; 

2.  The  adsorbed  atoms  encounter  each  other  on  the  solid, 

form  a product  molecule  and  emerge  into  the  gas 
phase. 

A key  assumption  in  this  research  is  that  the  energy  appearing  in 
vibration  of  the  product  molecules  (and  possibly  as  laser  radiation)  does  not 
come  from  the  solid  but  from  the  chemical  reaction  itself.  That  is,  as  a 
first  approximation  the  solid  is  treated  as  a rigid  substrate  that  alters 
the  energy  profile  from  gas  phase  reactants  to  gas  phase  products  but  does 
not  contribute  energy  to  the  reaction.  Energy  transfer  between  gases  and 
solids  is  often  a rather  inefficient  process  and,  we  speculate,  not  a likely 
source  of  sufficient  energy  to  lead  to  vibrational  excitation.  Therefore,  we 
look  to  the  reaction  itself  to  drive  the  vibrational  up-conversion,  and  in- 
clude the  dynamics  of  the  solid  only  as  a secondary  source  (or  sink)  for  the 
usable  energy.  Results  reported  in  Appendix  B confirm  that  the  energy  transfer 
between  the  reactants  and  the  solid  is  not  a major  factor  in  determining  the 
reaction  dynamics. 
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Using  this  assumption,  we  performed  an  extensive  series  of  calcu- 
lations for  the  recombination  step  (#2  above)  for  a variety  of  different 
systems.  This  was  performed  during  the  first  year  of  this  research  and  the 
details  are  re-printed  in  Appendix  A.  These  calculations  used  both  a rigid 
surface  and  a surface  permitting  energy  exchange.  As  described  in  Appendix  A, 
the  results  were  virtually  the  same;  showing  the  possibility  of  useful  popu- 
lation inversions  from  heterogeneously  catalyzed  reactions  if  certain  con- 
ditions (as  detailed  in  Appendix  A)  were  met. 

During  the  second  year  of  research,  the  adsorption  step  (#1)  was 
considered  and  two  effects  were  investigated.  First,  the  dynamics  of  adsorp- 
tion was  considered  in  order  to  estimate  the  amount  of  energy  transferred  to 
the  solid  (and,  hence,  lost  to  vibrational  up-conversion)  and  the  efficiency 
of  dissociative  chemisorption  for  a variety  of  potential  surfaces.  It  was 
found  that  dissociative  chemisorption  is  usually  a very  efficient  process  if 
thermodynamically  favorable.  That  is,  we  did  not  observe  serious  kinetic 
bottlenecks  to  dissociative  chemisorption  if  the  process  was  thermodynamically 
allowed  (e.g.,  exothermic  adsorption).  The  process  of  energy  transfer  to  the 
solid  was  more  troublesome.  All  studies  performed  to  date  assumed  the  surface 
was  rigid  (a  frozen  lattice  model)  or  the  surface  consisted  of  an  array  of 
uncoupled  enharmonic  oscillators  (an  Einstein  model).  We  observed  for  the 
Einstein  model  energy  transfers  to  the  solid  of  the  order  of  10%.  Therefore, 
it  became  important  to  make  this  estimate  more  accurate  and,  it  was  felt,  the 
use  of  uncoupled  oscillators  was  the  area  most  needing  improvement.  There- 
fore, the  second  phase  of  the  program  undertook  to  replace  uncoupled  oscillators 
with  generalized  Langevin  oscillators.  As  developed  by  Adelman  and  Doll^^^^ 
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the  generalized  Langevin  equations  (GLE)  are  a simple  extension  of  a model 
using  uncoupled  oscillators.  However,  rather  than  using  an  oscillator  with 
a conservative  force,  the  coupling  to  the  remainder  of  the  solid  is  intro- 
duced by  means  of  a damped  Brownian  oscillator.  That  is,  the  oscillator  is 
damped  to  approximate  energy  dissipation  into  the  solid  and  is  subject  to 
Brownian  random  forces  to  approximate  the  instantaneous  pumping  of  energy 
into  the  oscillator  from  the  random  motions  of  the  rest  of  the  solid.  The 
GLE  method  is  an  outgrowth  of  various  techniques  of  non-equilibrium  statis- 
tical mechanics  and  the  theory  of  stochastic  processes  as  developed  chiefly 
by  Kubo,^^^^  Mori,^^^^  and  Zwanzig.^^^^  The  GLE  method  has  proved  quite 
useful  for  gas-surface  collisions  and  is  likely  to  become  more  useful  in  a 
variety  of  areas  in  the  near  future. 

In  collaboration  with  Professor  S.  Adelman  of  Purdue  and  A.  Diebold 
(a  graduate  student  with  Adelman),  GLE  calculations  were  carried  out  for  com- 
parison with  Einstein  calculations  of  energy  transfer  to  the  solid.  The 
details  of  this  computation  and  the  results  are  reported  in  Appendix  C.  To 
summarize,  GLE  calculations  predict  increased  energy  loss  into  the  solid  when 
compared  to  Einstein  calculations,  by  roughly  a factor  of  2.  Therefore, 
perhaps  20%  of  the  reaction  energy  could  be  lost  to  the  solid  in  unfavorable 
circumstances  and  care  must  be  taken  to  minimize  this  effect. 


9 


IV. 


CUMULATIVE  LIST  OF  WRITTEN  PUBLICATIONS 


First  Year 

1.  "Pattern  Recognition  in  Chemical  Dynamics",  J.  H.  McCreery 
and  G.  Wolken,  Jr.,  Chem.  Phys.  Letters  469  (1977). 

2.  "A  Test  of  the  First  Order  Distorted  Wave  Born  Approximation 
in  Gas-Solid  Energy  Transfer",  Y.  W.  Lin,  S.  A.  Adelman, 

G.  Wolken,  Jr.,  Surface  Science  66,  376  (1977). 

3.  "Atomic  Recombination  Dynamics  on  Solid  Surfaces:  Effects  of 
Various  Potentials",  J.  H.  McCreery  and  G.  Wolken,  Jr.,  J. 
Chem.  Phys.  §1,  2551  (1977). 


Second  Year 

1.  "The  Energetics  of  Diatom/Solid  Dissociative  Adsorption", 

A.  C.  Diebold  and  G.  Wolken,  Jr.,  Surface  Science  (in  press). 

2.  "Vibrational  Excitation  From  Heterogeneous  Catalysis", 

G.  D.  Purvis  III,  M.  J.  Redmon,  and  G.  Wolken,  Jr.,  submitted 
to  J.  Phys.  Chem. 


V. 


LIST  OF  PROFESSIONAL  PERSONNEL  ASSOCIATED  WITH  THIS  RESEARCH 


Dr,  George  Wolken,  Jr.,  Associate  Section  Manager,  Chemical  Physics 
Group,  Principal  Investigator,  Battelle,  Columbus  Laboratories. 

Dr.  Jane  Hylton  McCreery,  Information  Systems  Section,  Battelle, 

Col umbus  Laboratori es . 

Dr.  George  Purvis  III,  Chemical  Physics  Group,  Battelle,  Columbus 
Laboratories. 

Dr.  Michael  J.  Redmon,  Chemical  Physics  Group,  Battelle,  Columbus 
Laboratories. 

Mr.  Alain  C.  Diebold,  Purdue  University  (graduate  student). 


n 


TRAVEL 


G.  Wolken,  Jr.,  discussion  leader,  Gordon  Conference  on  Atomic 
and  Molecular  Interactions,  July  25-28,  1978. 


12 


I 


REFERENCES 

1.  J.  C.  Polanyi,  Acc.  Chem.  Res.  5^,  161  (1972). 

2.  "Handbook  of  Chemical  Lasers",  Ed.  R.W.F.  Gross,  J.  F.  Bott  (John 
Wiley,  New  York,  1976),  Chapter  1. 

3.  K.  L.  Kompa,  "Chemical  Lasers"  (Springer-Verlag,  New  York,  1973). 

4.  Aviation  Week  and  Space  Technology:  August  21,  1978,  p.  38. 

5.  Aviation  Week  and  Space  Technology:  August  7,  1978,  p.  14. 

6.  J.  C.  Polanyi,  0.  Chem.  Phys.  84,  347  (1961). 

7.  J.V.V.  Kasper  and  G.  C.  Pimentel,  Phys.  Rev.  Lett.  1£,  352  (1965). 

8.  P.  J.  Kuntz  in  "Dynamics  of  Molecular  Collisions  - Part  B",  W.  H.  Miller, 
ed.  (Plenum  Press,  New  York,  1976),  Chapter  2. 

9.  B.  Halpern  and  D.  E.  Rosner,  J.  Chem.  Soc.  Faraday  Trans.  I,  74,  1883 
(1978). 

10.  D.  E.  Rosner,  private  communication. 

11.  S.  A.  Adelman  and  J.  D.  Doll,  J.  Chem.  Phys.  M,  2375  (1976). 

12.  R.  Kubo,  Rep.  Prog.  Theor.  Phys.  255  (1966). 

13.  H.  Mori,  Prof.  Theor.  Phys.  423  (1965). 

14.  R.  Zwanzig,  0.  Chem.  Phys.  1173  (1960). 


f 

I 

I 


I 

I 

I 


f - 

r 


j 


I 

\ 


11 


I 

I 

I 

I APPENDIX  A 

I "ATOMIC  RECOMBINATION  DYNAMICS  ON  SOLID  SURFACES: 

] EFFECT  OF  VARIOUS  POTENTIALS" 


I 


S' 

I 

1 


I 


I 

I 

I 

I 

I 

i 


I 

1. 

I 

1 

I 


1. 

I 


[ 


f 

I 

( 


Atomic  recombination  dynamics  on  solid  surfaces:  Effect  of 
various  potentials*^ 

Jane  Hylton  McCreery  and  George  Wolken,  Jr. 

Baitelle  Columbus  Laboratories,  Columbus.  Ohio  43201 
(Received  II  May  1977) 

* 

A model  potential  for  gas-solid  interactions  has  been  used  to  investigate  the  dynamics  of  recombination 
of  two  atoms  initially  adsorbed  on  a solid  surface.  In  the  spirit  of  Polanyi's  investigation  into  (he  eflect  of 
the  potential  energy  surface  on  the  dynamics  of  gas-phase  reactions,  a range  of  gas-solid  potential  energy 
surfaces  has  been  constructed.  Classical  trajectories  have  been  used  to  study  the  dynamics  of  reactions  on 
(hose  surfaces.  It  has  been  found  that  many  of  the  rules  postulated  by  Polanyi  for  energy  requirements 
and  disposal  mechanisms  for  gas-phase  systems  are  applicable  ilso  to  the  case  of  recombination  of 
adsorbed  atoms  to  form  a gas-phase  molecule.  Previous  work  assumed  a rigid  surface  providing  a static 
background  potential  in  which  the  adsorbed  atoms  moved.  An  extension  of  this  model  is  described  in  I 

which  the  rigid  surface  restriction  is  relaxed  and  one  or  more  surface  atoms  are  allowed  to  move  j 

interacting  with  the  adsorbed  atoms.  Using  this  potential  the  rigid  surface  model  is  shown  to  be  a good  i 

approximation  for  describing  many  aspects  of  recombination  dynamics.  j 


1.  INTRODUCTION 

A large  part  of  gas-phase  chemical  dynamics  is  con- 
cerned with  specific  energy  requirements  and  specific 
energy  disposal  mechanisms  in  chemical  reactions.  It 
is  now  abundantly  clear  that  not  all  types  of  energy  are 
equally  efficient  In  promoting  endothermic  reactions, 
nor  are  all  channels  of  energy  disposal  equally  efficient 
in  exothermic  reactions. ' To  assist  In  understanding  the 
microscopic  details  of  such  energy  specificity,  Polanyi 
and  co-workers  initiated  a series  of  model  calculations.  ‘ 
Much  of  their  work  used  several  simplifying  assump- 
tions: (a)  classical  trajectories  were  used;  (b)  atom- 
diatom  exchange  reactions  were  studied;  and  (c)  purely 
empirical  potentials  were  used  allowing  them  to  model 
a wide  range  of  realistic  <if  not  actually  existing  in  na- 
ture) potential  energy  surfaces.  The  resulting  “propen- 
sity rules”  have  found  wide  applicability  in  chemical  dy- 
namics. Most  of  the  essential  features  of  these  propen- 
sity rules  have  survived  the  test  of  more  refined  calcu- 
lations. The  great  utility  of  these  rules  lies  precisely 
in  their  qualitative  nature  and  broad  applicability  (e.  g. , 
attractive  potential  surfaces  lead  to  vibrational  exclta- 
tic  of  the  products).  Therefore,  it  is  not  surprising 
that  relatively  simple  calculations  proved  sufficient  to 
obtain  the  broadest  principles  of  energy  specificity. 

We  have  recently  developed  a model  potential  to  de- 
scribe the  interaction  of  diatomic’  and  trlatomic*  mole- 
cules with  rigid  (frozen-lattice)  solid  surfaces.  Using 
these  potentials,  classical  trajectories  have  been  calcu- 
lated to  study  selected  features  of  the  dynamics  of  atomic 
recombination,’  atom-adsorbed  atom  collisions,*  mole- 
cule-surface collisions,  ^ and  molecule-adsorbed  atom 
collisions.*  The  thrust  of  this  work  has  been  to  elucidate 
what  novel  dynamical  effects  (1.  e. , energy  specificity) 
may  he  present  in  chemical  reactions  occurring  on  solid 
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surfaces  and  to  what  extent  the  propensity  rules  devel- 
oped for  gas-phase  reactions  are  also  applicable  to  surface 
reactions.  It  Is  well  known  (and  of  considerable  econom- 
ic importance)  that  solid  surfaces  often  drastically  af- 
fect the  kinetics  of  chemical  reactions.  Although  reac- 
tion rates  are  considerably  affected,  relatively  little  at- 
tention has  been  paid  to  how  dynamics  (1.  e. , specific 
energy  requirements)  are  changed  by  the  use  of  hetero- 
geneous catalysts.  Computer  simulations,  in  the  spirit 
of  Polanyi's  work  In  gas-phase  dynamics,  are  well 
suited  to  exploring  this  question. 

For  the  Interaction  of  a diatomic  molecule  with  a solid, 
a first  approximation  could  be  to  use  gas-phase  informa- 
tion for  A4-(B,  C)  collisions,  but  permit  one  atom  to  be- 
come infinitely  massive,  thus  mimicking  a solid.  With 
this  assumption,  one  would  expect  heterogeneous  reac- 
tion dynamics  to  be  simply  a subset  of  gas-phase  dynam- 
ics. However,  there  are  two  complicating  factors:  (1) 
When  energy  transfer  to  or  from  the  solid  is  present, 
the  dynamics  could  be  affected  In  ways  not  predicted  by 
simple  A-f  (B,  C)  dynamics.  (This  degree  of  freedom 
was  not  present  in  our  previous  work  employing  the  rigid 
surface  model.  However,  surface  dynamics  is  included 
In  certain  of  the  calculations  described  below,)  (2)  A 
more  fundamental  difference,  present  also  in  the  rigid 
surface  approximation,  is  that  the  size  (as  well  as  the 
mass)  of  one  collision  partner  becomes  infinite.  There- 
fore, it  is  possible  (and  occurs  frequently  in  our  calcu- 
lations) that  two  atoms  can  both  be  close  to  the  solid,  but 
arbitrarily  far  from  each  other.  This  introduces  an  ad- 
ditional asymptotic  scattering  channel  which  has  no  good 
analogue  in  gas-phase  collisions  of  simple  molecules. 
(The  extent  to  which  uhimolecular  decomposition  or  col- 
lision complexes  of  large  molecules  mimic  gas-solid 
reaction  dynamics  is  not  clear.  Certainly,  gas-solid 
collisions  and  atom-diatom  collisions  are  two  limiting 
cases  for  which  interesting  dynamical  effects  could  be 
studied.  There  are  numerous  intermediate  cases.)  As 
discussed  below,  in  spite  of  these  complications,  many 
of  the  Polanyi  rules  for  simple  gas-phase  reactions  can 
be  formulated  in  a way  appropriate  for  the  class  of  sur- 
face reactions  studied  below. 


Tht  Journal  of  CheniMl  Phvtici.  Vol.  87.  No.  6.  IS  Seotimntf  1977  CoDvriaht  t)  1977  Ameticin  Initituw  of  Physio  2581 


r 

( 

'i 

P'' 


2552 


J.  H.  McCreery  and  G.  Wolken,  Jr.;  Atomic  recombination  on  surfaces 


In  the  present  studj’,  we  investigate  the  dynamics  of 
recombination  of  two  atoms  initially  adsorbed  on  a solid 
surface.  In  the  spirit  of  Polanyi’s  investigation  Into  the 
effect  of  the  potential  energy  surface  on  the  reaction  dy- 
namics for  gas-phase  systems,  ‘ we  use  a range  of  gas- 
sol  id  potential  energy  surfaces  generalizing  our  previous 
studies  of  recombination  dynamics.’  Also,  we  present 
preliminary  results  for  heterogeneous  reaction  dynamics 
in  which  the  surface  atoms  themselves  are  permitted  to 
move.  For  the  cases  studied,  no  significant  changes  re- 
sult from  motion  of  the  surface  atoms.  Section  11  gives 
a brief  overview  of  the  propensity  rules  developed  for 
gas-phase  collisions.  Section  ni  reviews  the  model  po- 
tential for  gas-surface  collisions  and  describes  the  range 
of  potentials  used  here.  Section  IV  presents  results  and 
discussion  for  the  rigid  surface  model.  The  generaliza- 
tion to  a moving  surface  is  given  In  Sec.  V,  and  prelimi- 
nary results  are  presented. 

II.  OVERVIEW  OF  GAS-PHASE  PROPENSITY  RULES 

For  exothermic  reactions  A+BC  — AB  + C with  exo- 
thermicin*  in  the  range  30-50  kcal  mole"'  and  activation 
barriers  typically  around  a few  kcal  mole*’,  It  was  shown 
that  attractive  potential  surfaces,  where  the  energy  Is 
released  as  A approaches  B,  lead  to  vibrational  excita- 
tion of  the  product  molecules;  repulsive  surfaces,  where 
the  energy'  Is  released  as  AB  separates  from  C,  give 
rise  to  rotation  and  translation  of  the  products.  An  ex- 
ception to  this  was  found  to  occur  for  repulsive  surfaces 
when  the  attacking  atom  A was  very  much  heavier  than 
atom  C.  This  leads  to  "mbted  energy  release”  dynamics 
and  product  vibrational  and  rotational  excitation.  Com- 
plex collisions,  where  more  than  one  encounter  between 
reagents  or  products  takes  place,  were  also  found  to  de- 
viate from  the  rule.  Such  secondary  encounters  mainly 
occur  on  attractive  surfaces  and  tend  to  reduce  the  prod- 
uct vibrational  excitation  on  highly  attractive  surfaces. 

Polanyi  and  co-workers  also  examined  the  effect  of  the 
position  of  the  crest  of  a barrier  for  thermoneutral  re- 
actions. It  was  found  that  “a  barrier  along  the  approach 
coordinate  is  most  efficiently  surmounted  by  motion 
along  the  approach  coordinate  (reagent  translation) 
whereas  a barrier  along  the  coordinate  of  separation  is 
most  efficiently  surmounted  by  motion  along  that  coor- 
dinate (reagent  vibration)."*  The  converse  was  found  to 
be  true  for  energy  hollows  rather  than  barriers.  These 
effects  were  found  to  be  independent  of  the  masses  of  the 
particles. 

For  endothermic  reactions  where  the  crest  of  the  bar- 
rier is  located  in  the  exit  valley,  it  was  found  that  vibra- 
tional energy  in  the  bond  under  attack  was  necessary  for 
successful  reaction.  An  exception  to  this  was  when  C had 
a much  s.-natler  mass  than  A or  B,  In  which  case  too 
much  vibration  In  the  reagents  markedly  reduced  the 
probability  of  reaction. 

These  results  were  concluded  from  series  of  classical 
trajector>’  calculations  on  London- Ey ring- Polanyl-Sato* 
(LEP?)-type  potential  functions  for  gas-phase  triatomlc 
systems.  The  aim  of  the  present  work  is  to  Investigate 
in  a similar  manner  the  effect  of  the  form  of  the  poten- 
tial energy-  function  on  the  dynamics  ot  reactions  on  solid 


surfaces.  Reactions  on  solid  surfaces  can  be  thought  of 
In  three  steps;  (1)  adsorption  of  the  gas  on  the  solid; 

(2)  possible  equilibrium,  or  partial  equilibration  of  the 
adsorbed  particles  with  the  solid;  (3)  recombination  and 
desorption  of  a gas  molecule.  If  the  adsorption  and  re- 
combination steps  are  fast  compared  to  the  adsorbate- 
lattice  relaxation  time.  Step  2 can  be  neglected.  The 
circumstances  under  which  this  occurs  are  the  circum- 
stances under  which  the  rigid  surface  model  should  be 
applicable.  This  is  discussed  more  fully  In  Sec.  V. 

The  purpose  of  the  present  study  is  to  attempt  to  mod- 
el the  Internal  vibration-rotation  states  of  diatomic 
molecules  as  they  are  formed  and  desorb  from  the  sur- 
face and  to  understand  which  features  in  the  potential 
lead  to  what  internal  states.  As  Polanyi  first  discussed 
for  gas-phase  exothermic  reactions, ' circumstances  ex- 
ist in  which  useful  amounts  of  the  reaction  energy  appear 
as  vibration  in  the  products.  A further  goal  of  the  pres- 
ent computer  simulations  is  to  explore  (admittedly,  in  a 
very  crude  way)  possible  situations  under  which  similar- 
ly useful  vibrational  populations  might  be  obtained  from 
heterogeneously  catalyzed  reactions.  The  function  of  a 
Catalyst  is  to  change  the  energetics  of  the  reaction  path 
to  accelerate  the  formation  of  products.  Similarly,  a 
change  in  the  energetics  of  the  reaction  path  will  affect 
the  distribution  of  the  available  energy  in  the  products. 

It  is  this  modified  energy  distribution  we  will  attempt  to 
model. 

We  envision  attempting  to  simulate  the  following  pro- 
cess; consider  a (hypothetical)  reaction  of  gas  molecules 
that  is  highly  exothermic  but,  for  whatever  reason,  does 
not  produce  the  desired  state  distribution  of  products. 

A catalyst  Is  used  that  alters  the  energetics  of  the  reac- 
tion path  and  produces  some  other  product  state  distri- 
bution. The  available  energy  comes  from  the  gas-phase 
reagents  with  the  solid  remaining  chemically  and  dynam- 
ically Inert  throughout  the  reaction.  (We  will  relax  the 
restriction  that  the  solid  be  strictly  rigid  in  Sec.  V,  but 
we  will  assume  that  virtually  all  the  available  energy 
comes  from  the  gas-phase  reagents.  The  extent  that  the 
solid  is  an  energy  sink  under  these  conditions  will  be 
discussed  In  Sec.  V.)  While  the  temperature  of  the  solid 
Can  often  affect  its  catalytic  activity,  it  is  usually  as- 
sumed that  the  primary  mechanism  is  that  of  a hot  sur- 
face effectively  lowering  activation  barriers  more  than  a 
cold  surface.  Since  this  does  not  require  gas- surface 
energy  transfer.  It  too  can  be  treated  with  either  the 
rigid  surface  or  the  moving  surface  model. 

Hence,  we  are  concerned  with  basically  a two-step 
process:  adsorption  followed  by  recombination  and  de- 
sorption with  a net  energy  release  into  the  products  and 
possibly  into  the  solid.  The  present  work  considers  only 
the  recombination  processes  to  assess  how  the  potential 
surface  affects  the  product  state  distributions.  The  dy- 
namics of  the  adsorption  process  determine  the  Initial 
conditions  for  the  recombination  (i.  e. , the  initial  mo- 
menta of  the  adsorbed  atoms  when  they  collide)  and 
therefore  is  one  step  further  removed  from  the  formation 
of  products.  The  effect  of  various  initial  momenta  on 
product  state  distributions  Is  the  subject  of  a future 
study. 
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FIG.  1.  Geometry  of  the  W(OOl)  surface  indicating  the  singly 
coordin.-ited  (1C>V)  adsorption  site  and  the  2CS,'  oCS 
The  origin  of  coordinates  is  at  a ICA’  site  with  the  x-y  plane 
parallel  to  the  surface  and  the  x axla  directed  towards  a neigh- 
bori.ng  IC.V  site.  To  study  recombination  dynamics,  two  atoms 
are  given  random  initial  momenta,  subject  to  the  restriction 
that  the  (x,y)  components  of  momenta  lie  within  the  square  in- 
dicated. This  neglects  “catch-up"  coUlalona.  a«5.97a. u. 

We  have  also  considered  the  question:  how  many,  it 
any,  of  the  "Polanyi  rules”  for  gas-phase  dynamics  are 
applicable  also  to  gas-surface  dynamics.  In  particular 
to  recombination  dynamics?  For  example,  in  an  exo- 
thermic recombination  reaction  is  it  true  that  an  attrac- 
tive potential  surface  will  give  rise  to  vibrational  excita- 
tion of  the  product  molecule?  The  potential  functions 
used  throughout  this  work  have  the  form  of  the  modified 
LEPS  potential  that  has  been  used  in  previous  studies  of 
gas-solid  interactions.’'^  By  varying  the  parameters  in 
the  potential  function  we  generate  a wide  variety  of  po- 
tential surfaces.  We  use  both  equal  mass  and  unequal 
mass  combinations.  We  have  relaxed  the  rigid  surface 
approximation  in  part  allowing  some  of  the  surface  atoms 
to  move,  and  we  examine  the  effect  of  the  surface  motion 
on  the  recombination  dynamics. 

III.  INTERACTION  POTENTIALS;  RIGID  SURFACE 

The  potential  functions  used  were  generated  from  the 
diatom-solid  surface  modified  LEPS  potential  that  has 
been  used  in  previous  work.  This  has  the  form 

I'ttp*  * £/i  + + I/j  -[A  \ + (Aj+Aj)’  - Ai(i4j+Aj) 

(1) 

where 

{<3 + A()  exp[-  2a, (y,  - r,o) ] 

-(2  + 6A,)exp[-a,(r,-r,J]},  (2) 


^ * 3a,)  expl-  2a, (r,  - rjo)] 

-(6+2A,)exp[-a,(r, -r,o)]}  , (3) 

and  Di,  a,,  and  rio  are  the  dissociation  energy,  Morse 
parameter,  and  equilibrium  distance  for  the  ith  two-body 
interaction.  We  take  i = 1 to  be  the  atom-atom  interac- 
tion and  assume  a typical  Morse  curve  to  obtain  the 
parameters  f),,  a,  and  y^.  A,  is  used  as  a parameter  to 
generate  different  potential  surfaces,  i > 2,  3 correspond 
to  the  atom- surface  interactions,  and  it  is  assumed  that 
these  are  the  same  for  both  atoms.  To  account  for  the 
structure  of  the  solid  in  the  atom-surface  interactions 
we  require  the  parameters  D,  a and  yg  to  be  functions  of 
X and  y,  where  the  x-y  plane  is  parallel  to  the  plane  of 
the  solid  surface.  The  symmetry  assumed  for  the  sur- 
face is  that  of  the  (100)  face  of  a bcc  solid  and  the  geo- 
metry is  that  of  tungsten  (see  Fig.  1).  Following  the 
previous  work  of  Hj  interacting  with  the  (100)  face  of 
tungsten,  we  assume  the  existence  of  three  possible 
binding  sites,  the  on-top  (ICN),  bridge  {2CS),  and  hole 
(SCAT)  sites.  We  take 

f>(x,y)  =f7o[l+aQ(x,y)],  (4) 

To(x,  y)  =Z,g[l  + cP(x,  y)],  (5) 

where 

Q(.x,  y)  **|^coo^^^j + cos^-^^  j 

-a[cOs(^)-i][cos(^)-i],  (6) 


P(x,y).*[cos(^)^cos(^)] 


-b[cos(^)-i][cos(^)-i]. 

a is  the  tungsten  nearest  neighbor  distance  (>  5. 97  a,  u.). 
k is  chosen  to  be  either  0 or  1;  in  the  previous  work  it 
was  always  assumed  to  be  1.  The  effect  of  - 0 is  to 
Introduce  channels  into  the  egg-:a;U>n  shape  of  the  sur- 
face by  making  the  on-top  and  bridge  sites  Identical  for 
hydrogen  binding.  The  Morse  parameter  a is  chosen  as 
before  to  be  [0. 02894/f)(x,  y)]*'*,  where  D is  in  atomic 
units  of  energy  giving  a in  atomic  units  of  inverse  length 


We  have  a total  of  nine  parameters  with  which  to  gen- 
erate different  potential  surfaces,  namely,  Z)g,  5,  k.  A, 
Z„,  e,  B,  A„  A(«Aj«Aj).  The  parameters  were  ad- 
justed to  give  nine  different  potential  energy  surfaces 
with  varying  barriers  and  modes  of  energy  release. 
These  parameters  are  given  in  Table  I.  Equipotential 


TABLE  1.  Parameters  used  in  generating  the  varloui  potentials.- 


Binding  energy  (eV) 


Potential 

c 

Zm 

6 

D, 

* 

A 

B 

A, 

Aj* 

ICN 

ZCN 

5CA’ 

I 

1.0 

3.12 

1.0 

0. 1102 

0.0 

0. 1067 

0.2051 

0. 147 

0.03 

3.0 

3.0 

1.72 

U 

0.2112 

2.19 

0. 0474 

0. 1007 

1.0 

1.463 

0.653 

0.147 

0.  03 

3.0 

2.74 

1.72 

lU 

0.2112 

2.19 

0.0474 

0. 1007 

1.0 

0. 1463 

0.653 

-0. 1 

-0.  1 

3.0 

2.74 

2.40 

IV 

1.0 

3.12 

1.0 

0.0S51 

0.0 

0.1067 

0.2051 

0.147 

0.03 

1.5 

1.5 

0.  86 

V 

1.0 

3.12 

1.0 

O.OSSl 

0.0 

0. 1067 

0.2031 

0.3 

0.3 

1.5 

1.5 

0.  86 

VI 

0.2112 

2. 19 

0.  0474 

0.  0503 

1.0 

1.463 

0.653 

0.  147 

0.03 

1.5 

1.37 

0.  86 

vn 

1.0 

3.12 

1.0 

0.0551 

0.0 

0. 1067 

0.2051 

0.0 

-0.3 

1.5 

1.5 

0.  SC 

vm 

1.0 

3.  12 

1.0 

0.0551 

0.0 

0. 1067 

0.2051 

0.  147 

0.3 

1.5 

1.5 

0.  86 

L\ 

0.2112 

2. 19 

-0.  1S81 

0. 0503 

1.0 

-0.3725 

0.653 

0.147 

0.03 

0.  SC 

1.37 

1.3 
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FIG.  2.  Equipotcntial  contour*  (in  eV)  (or  th*  approach  o(  H;  to  W(OOl)  (or  the  di((erent  potential  turtace*  itudied.  The  contour 
level*  are  equally  epaced  (or  any  given  potential  *ur(ace.  The  H-H  bond  1*  parallel  to  the  plane  at  the  euifaee.  The  center  o( 
mae*  o(  the  molecule  1*  perpendicularly  above  a 2CN  alta.  X Indicate*  the  poeltlon  o(  the  reactant*,  two  atom*  ad*orb*d  *t  ad- 
jacent IC.V  site*  {SCN  (or  IX).  • locate*  our  choice  of  the  “traneitlon  eite.” 


contour  plots  are  given  in  Fig,  2.  These  are  (or  a hy-  surfaces,  it  Is  stretching  along  the  y axis.  The  dashed 

drogcn  molecule  approaching  the  surface  with  the  axis  line  Indicates  the  minimum  energy  reaction  path  from 

parallel  to  the  plane  of  the  surface.  The  midpoint  of  the  the  reactants , two  atoms  adsorbed  on  the  surface,  to 

bond  is  perpendicularly  approaching  a bridge  (or  2CN)  products,  a molecule  far  from  the  surface.  The  X 

site  and  the  bond  is  stretching  symmetrically  In  the  dl-  marks  the  initial  position  of  the  atoms  In  the  ICS  or 
reetton  of  the  atomic  adsorption  site  of  greatest  stabll-  iCS  sites.  The  energy  along  the  reaction  path  Is  shown 

ity.  This  is  the  on-top  (IC.V)  site  for  all  the  potential  in  Pig,  3.  The  dashed  line  parallel  to  the  ordinate  indl- 

.surf.-icbs  c-xccpt  DC,  where  the  hole  site  (iCS)  is  the  cates  the  position  along  the  reaction  path  of  an  arbitrari- 

most  stable  site.  With  the  coordinate  system  of  Fig.  1,  ly  defined  "transition  state"  between  reacUnts  and  prod- 

the  bond  in  potential  DC  is  stretching  parallel  to  the  x uett.  While  this  transition  state  is  not  well  defined  for 

axis  along  the  line  y >a/2,  while  for  the  other  potential  these  molecule- surface  potential  functions,  it  is  a con- 


J.  H.  McCreery  and  G.  Wolken,  Jr.:  Atomic  recombination  on  surfaces 


2555 


FIG.  3.  Energy  profiles  along  the  reaction  paths  for  the  potential  surfaces  given  In  Fig.  2.  The  dashed  line  marks  the  region  of 
'rtnsltlon  from  reagents  to  products  (denoted  by  • In  Fig.  2). 


cept  useful  (or  defining  the  boundary  between  entrance 
and  exit  channels  and  therefore  plays  a role  in  any  dis- 
cussion of  the  effects  of  barrier  position  in  relation  to 
these  channels.  Consequently,  we  have  chosen  to  deter- 
mine the  transition  state  region  approximately  for  each 
potential  surface.  In  potentials  IV,  V,  and  VIII  the  dis- 
tance along  the  reaction  path  of  the  transition  state  re- 
gion Is  easily  defined  since  the  reaction  path  makes  quite 
a sharp  turn  at  that  point.  For  the  other  potentials,  the 
exact  position  is  not  clearly  defined  but  we  chose  it  to  be 
approximately  half-way  around  the  curve  in  the  reaction 
path  in  the  transition  state  region. 

Potentials  I-m  are  endothermic  with  endothermictty 
of  approximately  1, 26  eV.  We  assume  that  the  overall 
reaction  catalyzed  by  the  solid  surface  is  exothermic. 
However,  it  may  occur  that  the  adsorption  step  releases 
more  energy  than  the  total  reaction.  Therefore,  the  re- 
combination will  occur  between  atoms  with  considerable 


translational  energy,  but  on  an  uphill  reaction  path.  Po- 
tentials I-in  are  intended  to  model  this  possibility.  Po- 
tential I has  the  crest  of  the  barrier  (of  height  1. 39  eV) 
in  the  exit  channel.  The  on-top  and  bridge  sites  are 
identical  for  this  case  and  are  more  stable  than  the  hole 
site,  n also  has  the  crest  of  the  barrier  (of  height  1. 39 
eV)  in  the  exit  channel.  The  adsorption  sites,  in  order 
of  decreasing  stability,  are  on-top,  bridge,  and  hole. 

For  gas-phase  reactions  on  potential  surfaces  of  this 
form,  vibrational  energy  in  the  bond  under  attack  is  nec- 
essar>’  for  successful  reaction.  It  is  not  immediately 
obvious  what  analogous  condition  on  the  gas- solid  poten- 
tials will  be  necessary  for  successful  recombination. 
Potential  in  is  also  ondothermiv.  but  has  the  crest  of  the 
barrier  (of  height  1.99  oVi  at  more  or  less  the  transition 
state  position. 

The  remaining  potential  functions  are  exothermic  with 
exothermlclty  of  approximately  1.73  eV.  Potential  IV 
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has  no  barrier  to  molecule  formation  and  desorption, 
and  the  energj"  release  occurs  about  equally  in  both  en- 
terance  and  exit  channels.  For  gas-phase  reactions, 
this  mixed  energy  release  would  lead  to  vibration,  rota- 
tion, and  translation  of  the  products.  Potential  IV  also 
has  no  barrier  from  reagents  to  products  but  for  this 
case  much  of  the  energy  release  takes  place  in  the  en- 
trance channel.  This  is  an  attractive  surface  and  for 
gas-phase  systems  would  lead  to  vibrationally  excited 
product  ;nolecules.  Potential  VI  has  a small  barrier  of 
height  0. 10  eV  entirely  in  the  entrance  channel.  As  with 
rv  about  half  the  energy  release  occurs  in  each  channel. 
If  an  analogy  can  be  drawn  with  the  gas-phase  reactions, 
the  barrier  will  be  most  efficiently  surmounted  by  trans- 
lational energy  , VII  also  has  a barrier  (of  height  0. 33 
eV)  with  crest  just  in  the  entrance  channel.  The  energy 
release  occurs  mostly  in  the  exit  channel.  This  Is  a re- 
pulsive surface  requiring  vibrational  energy  for  success- 
ful reaction  in  the  gas-phase  case.  VUI  has  no  barrier 
and,  as  with  V,  has  nearly  all  the  energy  released  in  the 
entrance  channel.  Potential  DC  is  generated  from  the 
same  data  as  VI  except  that  the  relative  stability  of  the 
adsorption  sites  is  reversed,  the  hole  site  being  the 
most  stable  followed  by  the  bridge  and  then  the  on-top 
sites.  This  also  leads  to  a mixed  energy  release  type  of 
surface  when  the  two  atoms  are  initially  positioned  at 
adjacent  hole  sites. 

IV.  CALCULATIONS 

For  each  potential  surface  we  computed  a series  of 
classical  trajectories  designed  to  simulate  the  recom- 
bination of  tivo  adsorbed  atoms.  Initial  conditions  for 
those  trajectories  have  been  described  in  previous  work 
on  recombination  dynamics.  ® The  atoms  were  positioned 
initially  at  adjacent  IC.V  sites  (5C.V  sites  for  potential 
DO.  The  total  initial  kinetic  energy  was  fbced  but  was 
randomly  distributed  between  the  two  atoms.  The  initial 
momenta  were  chosen  w’ith  random  directions  but  re- 
stricted so  the  r and  v components  were  directed  within 
the  square  shown  in  Fig.  1.  The  total  initial  kinetic  en- 
ergy for  the  system  was  chosen  to  be  3, 50  eV  for  the  en- 
dothermic potential  surface  and  0.  50  eV  for  the  exother- 
mic surfaces.  These  values  are  sufficiently  large  to 


TABLE  11.  Mean  distribution  of  energy  in  product  molecules. 


Potential 

' Molecule 
formation 

c>  r 

C F 

* 

U)  Equal  masses 

IV 

33 

43 

IS 

42 

3o 

57 

23 

20 

VI 

33 

50 

10 

40 

VIII 

29 

61 

20 

19 

IX 

32 

49 

5 

46 

(b)  Lisht- 

-heav>’  mass  combination 

IV 

24 

46 

12 

42 

39 

44 

34 

22 

VI 

37 

52 

13 

35 

vi:i 

2(i 

4«/ 

32 

23 

fX 

4ft 

30 

20 

50 

TABLE  111.  Product  vibrational  slate  dis- 


tribution summed  over  rotational  sublevels. 


Potential  t’«0 

1 

2 

3 

4 

5 

(a)  Equal  masses 

IV 

IS 

39 

30 

12 

1 

0 

11 

17 

36 

34 

3 

0 

VI 

S 

36 

39 

15 

2 

0 

VIII 

1 

20 

48 

26 

5 

0 

L\ 

14 

29 

32 

24 

1 

0 

(b)  Light- 

•heavy 

mass 

combination 

IV 

11 

20 

20 

30 

15 

4 

V 

15 

21 

26 

22 

7 

9 

VI 

10 

10 

23 

23 

23 

10 

vai 

IS 

23 

13 

13 

27 

8 

IX 

39 

OO 

14 

20 

3 

2 

surmount  any  barriers  and  to  allow  a wide  range  of  prod- 
uct states.  For  each  potential  surface  we  ran  trajec- 
tories for  atoms  of  equal  mass  (hydrogen)  and  for  one 
light  (mass  Ij  and  one  heavy  (mass  80)  atom.  Three 
hundred  trajectories  were  computed  for  each  case. 

For  the  endothermic  potential  surfaces  I-III,  very  few 
molecules  W'ere  formed  for  either  the  equal  mass  or  the 
light-heavy  mass  cases.  Thus,  without  running  a very 
large  number  of  trajectories  (as  was  done  in  Refs.  6 and 
7),  we  have  insufficient  numbers  of  product  molecules 
to  draw  any  statistically  meaningful  conclusions  about 
properties  of  the  product  molecules.  For  successful 
gas-phase  reactions  on  endothermic  potential  surfaces, 
vibrational  energy  In  the  bond  under  attack  is  necessary. 
There  is  no  direct  analogue  of  this  vibrational  energy  in 
the  gas-solid  recombination  dynamics,  but  motion  on  the 
solid  surface  does  not  appear  to  be  sufficient  for  success- 
ful recombination  on  such  potentials.  Of  course,  even 
a low  probability  of  reaction  per  collision  may  lead  to 
molecule  formation  and  desorption  after  a sufficient 
number  of  collisions.  This  would  require  a much  longer 
residence  time  on  the  surface  and,  hence,  a greater  pos- 
sibility for  gas— solid  energy  transfer.  Such  equilibra- 
tion with  the  surface  is  inconsistent  with  the  constraints 
of  our  model  and,  therefore,  actually  running  targe 
numbers  of  trajectories  is  not  likely  to  contribute  to  an 
understanding  of  the  processes.  We  conclude  that  endo- 
thermic desorption  processes  appear  to  be  relatively  in- 
efficient at  forming  molecules,  even  with  considerable 
kinetic  energy.  Increased  energy  transfer  with  the  solid 
is  expected  and  hence  increased  “leakage"  of  reaction 
energy  on  the  way  from  reactants  to  products  is  also  ex- 
pected. Therefore,  reactions  with  very  strongly  bound 
intermediates  are  not  likely  candidates  to  channel  energy 
efficiently  into  the  product  molecules. 

Tables  n and  III  give  the  results  of  the  trajectory  cal- 
culations on  the  exothermic  surfaces.  Potential  VH  is 
not  Included  since  in  this  case  also  an  Insufficient  num- 
ber of  molecules  were  formed.  This  Is  a repulsive  sur- 
face, having  most  of  the  energ>'  released  in  the  exit 
channel,  and  in  the  gas  phase  such  a surface  al.so  re- 
quires vibrational  energy*  in  the  bond  under  attack  for 
successful  reaction.  While  there  seems  to  be  no  gas- 
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solid  analoiiue  to  such  vibrational  energj',  we  have  again 
lound  that  translation  over  the  surface  (including  compo- 
nents of  motion  perpendicular  to  the  surface)  is  not  suf- 
ficient despite  the  exothermicity  of  the  reaction.  For  the 
remaining  exothermic  potential  surfaces  approximately 
onc-thu'd  of  the  trajectories  led  to  reaction  and  the  for- 
mation of  a product  molecule.  Tlie  total  energy  avail- 
able to  these  molecules  for  internal  excitation  (inciuding 
zero-point  energy!  and  translation  is  about  2.25  eV.  Ta- 
ble II  gives  the  mean  distribution  of  this  energy  between 
vibration,  rotation,  and  translation.  Table  ni  gives  the 
vibrational  state  distributions  summed  over  the  rota- 
tional sublevels  of  the  desorbed  product  molecules.  For 
the  case  of  hvo  hydrogen  atoms,  the  attractive  potential 
surfaces  V and  Vin  lead  to  considerable  vibrational  ex- 
citation, as  is  the  case  for  attractive  surfaces  in  gas- 
phase  reactions.  These  surfaces  also  lead  to  the  most 
rotational  excitation  of  the  products  (and  hence  the  least 
translatiOiial  energy,  consistent  with  the  attractive  po- 
tential surfaces).  Potential  rv"  is  a mbced  energy  re- 
lease surface  and  leads  to  just  about  equal  amounts  of 
energy  in  vibration  and  translation.  This  can  be  re- 
garded as  analogous  to  the  gas-phase  results  for  mixed 
energy  release  surfaces.  Potentials  VI  and  DC  are  also 
mtxed  energy  release  but  each  has  a small  barrier  in 
the  entrance  channel.  The  initial  translational  energy 
of  the  adsorbed  atoms  is  seen  to  be  adequate  to  sur- 
mount these  barriers,  again  consistent  with  gas-phase 
results.  VT  and  DC  differ  only  in  the  relative  ordering 
of  the  adsorption  sites  and  this  difference  seems  to  have 
little  effect  on  the  recombination  properties  of  the 
surfaces. 

The  light-heat7  mass  combination  has  the  effect  of 
smoothing  out  the  sharper  peaks  In  the  vibrational  state 
distribution.  For  the  mixed  energj'  release  potentials 
IV  and  VI  (which  for  H,  peak  around  (■«  1,  r«2)  this 
causes  the  mean  vibrational  energy'  to  Increase  very 
slightly.  For  the  other  potentials  which  peak  around  v 
>2,  r - 3,  the  mean  vibrational  energy  is  decreased  con- 
siderably by  this  smoothing.  The  proportion  of  energy 
that  occurs  as  translation  of  the  product  molecules  var- 
ies little  with  the  mass  combination.  For  the  llght- 
hoatw  case  this  means  a much  smaller  velocity  of  the 
molecule,  and  indeed  these  trajectories  took  much 
longer  to  compute  than  did  the  equal  (light)  mass  case. 
Since  the  translational  energy  is  roughly  the  same  for  the 
nvo  mass  cases,  the  change  in  mean  vibrational  energy 
is  reflected  in  a change  in  mean  rotational  energy.  The 
unequal  mass  case  gives  rise  to  increased  rotational  en- 
ergy for  the  attractive  potentials  surfaces  V and  VIII  and 
also  for  the  mixed  energy  release  surface  DC.  Poten- 
tials VI  and  DC  differ  only  in  the  relative  ordering  of  the 
adsorption  sites  and  give  similar  results  for  equal 
masses,  but  for  light- heavy  mass  combinations  DC  gives 
rise  to  significantly  less  vibrational  excitation  than  VI 
and  correspondingly  more  rotation  and  translation. 

V.  MOVING  SURFACE 

The  rigid  surface  model  was  previously  justified’  by 
the  aiiijn.ption  that  the  adsorlK‘d  atoms  encounter  each 
other  and  recombine  (or  not)  on  a time  scale  short  com- 


pared to  adsorbate-lattice  energy  transfer  times.  For 
light  adsorbates  and  hea\7  surface  atoms,  it  was  felt 
that  the  mass  difference  would  ensure  only  modest  ener- 
gy transfer.  Our  previous  estimates  were  2^  energy 
transfer  for  H atom-tungsten  interactions,'"  or  about  4(7 
for  Hj- tungsten  interactions.  However,  since  we  are 
using  heavier  adsorbate  atoms  in  some  of  the  present 
calculations,  it  is  appropriate  to  investigate  more  care- 
fully the  validity  of  the  rigid  surface  model.  Therefore, 
we  have  run  trajectories  with  the  rigid  surface  restric- 
tion relaxed  and  some  of  the  surface  atoms  allowed  to 
move. 

It  is  first  necessary  to  generalize  Vleps  (the  LEPS 
potential  for  the  interaction  of  a gas  molecule  with  a 
rigid  solid  surface).  In  addition  to  a gas-rigid  surface 
Interaction,  can  also  be  considered  as  the  inter- 

action of  a gas  molecule  with  the  lattice  sites  of  the 
solid.  If  we  allow  motion  of  the  solid  atoms,  at  any  in- 
stant those  atoms  may  or  may  not  be  at  a lattice  site. 
Therefore,  we  have  modified  Vleps  by  Introducing  cor- 
rection terms: 

(1)  J'p,  the  restoring  force  on  each  atom  of  the  solid, 
tending  to  return  it  to  its  lattice  site; 

(2)  I coRBj  to  account  for  the  change  in  the  gas- sur- 
face interaction  when  the  solid  atom  is  displaced  from 
its  lattice  site. 

We  have 

Vtpt«i  = Vleps+ V;,  + VjQpn  . (8) 

V;,  and  VcoRR  must  vanish  when  all  the  solid  atoms  oc- 
cupy their  lattice  positions.  We  have  used  pair  poten- 
tials for  Vj,  and  Vcork  t®  correct  Flips  approxi- 
mate way  to  allow  for  motion  of  the  solid  atoms.  As- 
suming that  Vp  binds  the  solid  atom  to  its  lattice  site 
constitutes  an  Einstein  model  of  the  solid.  Certainly, 
pairwise  atom-atom  interactions  can  be  used  for  Vr, 
and  are  commonly  used  to  simulate  properties  of  bulk 
solids, However,  we  feel  that  the  simplicity  of  the 
Einstein  model,  and  the  corresponding  savings  of  com- 
puter time,  is  justified  for  the  present  rough  estimates 
of  energy  transfer.  More  refined  energy  transfer  cal- 
culations, using  a generalized  Langevin  oscillator  model 
for  the  solid,  demonstrate  that  the  Einstein  model  is  a 
reasonable  first  approximation  for  short  collision  times. " 

The  model  does  not  include  the  nonadditive  correc- 
tions for  the  solid- sol  id  interactions,  or  the  nonadditive 
corrections  to  Vconif  the  restoring  force  Vp  we 
hat’e  used  a harmonic  potential  binding  an  atom  of  the 
solid  to  its  (fixed)  lattice  site.  To  account  for  the  change 
in  the  gas-solid  potential  due  to  displacement  of  the  sur- 
face atom  from  its  lattice  site,  wc  have  used  a pair  po- 
tential for  Vcopp  connecting  each  gas  atom  with  each 
surface  atom.  Vlcpp  already  contains  the  Interaction  of 
each  gas  atom  with  each  lattice  site,  so  to  avoid  includ- 
ing it  twice,  this  lattice  site  interaction  must  be  sub- 
tracted from  fcopp.  We  take 

Fcorr'  El!  hl-(/U,) 1.  '9' 

i ' 

where  .V,  is  ttie  number  of  gas  atoms,  .V,  is  the  number 
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TABLE  IV.  Mean  distribution  of  energy  in  product  molecules: 
(a)  rigid  surface;  (b)  two  surface  atoms  free  to  move. 


Potential 

% Molecule 
formation 

% 

% Energy 
transfer 

(a)  Equal  maasea 

IV 

(a) 

41 

41 

15 

44 

lb) 

72 

49 

11 

34 

6.32 

V 

ta) 

36 

60 

20 

20 

lb) 

78 

52 

26 

IS 

6.69 

VI 

(a) 

42 

51 

10 

39 

lb) 

65 

57 

8 

30 

5.24 

VIII 

u) 

30 

68 

18 

14 

lb) 

68 

62 

18 

13 

7.04 

(b)  Light-heavy  maae  combination 

IV 

la) 

21 

47 

13 

40 

lb) 

73 

51 

10 

34 

4.79 

V 

(a) 

40 

49 

30 

21 

lb) 

79 

37 

40 

19 

3.86 

VI 

<a) 

40 

51 

13 

36 

lb) 

69 

47 

17 

32 

5.89 

vm 

(a) 

26 

37 

40 

23 

lb) 

70 

46 

27 

23 

4.26 

of  solid  atoms  that  are  free  to  move,  R^,  is  the  distance 
between  the  gas  atom  and  the  solid  atom,  and  is  the 
distance  between  the  gas  atom  and  the  lattice  site  of  the 
solid  atom.  As  required,  if  the  surface  atom  occupies 
its  lattice  site,  A,.,  > A,.,  and  Vcomi  > 0.  The  total  gas- 
solid  potential  for  the  moving  surface  case  is 

• < • 

where  Vf  is  the  harmonic  restoring  force  for  solid  atoms. 

To  summarize,  the  revised  model  includes 

(1)  all  forces,  pairwise  and  nonpairwise,  to  describe 
the  interaction  of  a gas  molecule  with  a surface  when  the 
surface  atoms  are  fbced  at  their  lattice  sites; 

121  pairwise  corrections  for  the  motion  of  the  solid 
atoms  away  from  their  lattice  sites; 

(3)  pairwise  corrections  to  account  for  the  change  in 
the  gas-solid  potential  due  to  displacement  of  the  surface 
atoms  from  their  lattice  sites. 

We  have  used  this  potential  to  examine  the  validity  of 
the  rigid  surface  approximation  for  the  recombination 
studies.  We  require  the  harmonic  restoring  force  for 
the  solid  atoms,  and  the  gas-solid  atom  Morse  potential 
[ll'  in  Eq.  (9)].  For  the  harmonic  restoring  force  we 
need  an  estimate  of  the  Einstein  temperature  of  the  sur* 
face.  The  Einstein  temperature  can  be  estimated  from 
the  Delr.-e  temperature  either  by  requiring  the  best  sln- 
gie-oseillator  approximation  to  the  Debye  mode  density” 
or  by  e.s'panding  the  high  temperature  Debye  and  Einstein 
heat  capacities  and  matching  terms  through  T**.  In  both 
cases,  one  finds  the  Debye  temperature  to  be  smaller 
lli..;i  liif  EiaMvin  temperature  by  a factor  of  % 0. 0.  Also, 
in  Ol  der  to  correct  approximately  for  the  smaller  Debye 


temperature  of  the  surface  compared  to  the  bulk,  an  ex- 
tra factor  of  (2/3)  is  used.  Hence, 

0Erj'i¥.V  = (2/3)'/O'  (U) 

Using  = 400  "K  yields  8ET!i“iV  =•  210 °K.  This  gives 
a harmonic  force  constant  for  a tungsten  solid  atom  of 
0. 074  a.  u.  The  solid  atoms  were  assumed  initially  to 
be  at  their  equilibrium  positions  and  those  that  were  free 
to  move  were  given  an  initial  kinetic  energy  (the  same 
for  all  the  atoms)  equal  to  the  average  energy  of  the  har- 
monic oscillator  at  a surface  temperature  of  300  °K  with 
an  Einstein  temperature  of  210°K.  Then  each  solid  atom 
initially  has  kinetic  energy  of  0. 001B7  a.  u. ; the  initial 
direction  of  velocity  of  the  atoms  was  randomly  chosen 
so  that  there  was  equal  probability  of  it  being  within  any 
differential  solid  angle.  The  gas— solid  Morse  potential 
used  in  Fcokr>  *^(^-«)>  chosen  to  have  a>0.5123 
a.u.;  D^O.OSSa.u.  and  ro^  3. 12 a. u.  These  parame- 
ters correspond  to  the  hydrogen  atom-tungsten  ICV 
Morse  potential  for  the  surfaces  used. 

To  examine  the  validity  of  the  rigid  surface  model  In 
recombination  dynamics,  100  trajectories  were  computed 
for  each  of  potentials  IV,  V,  VI,  and  vni  and  for  both 
the  equal  (light)  mass  and  the  light-heavy  mass  combina- 
tion. The  Initial  conditions  of  the  adsorbed  gas  atoms 
were  chosen  to  be  the  same  In  these  trajectories  as  in 
the  first  100  trajectories  for  those  potentials  in  the  rigid 
surface  model.  Two  surface  atoms  were  allowed  to 
move,  those  being  the  two  atoms  in  the  ICV  sites  direct- 
ly below  the  adsorbed  gas  atoms  In  their  initial  configu- 
ration. The  Initial  directions  of  velocity  of  these  two 
surface  atoms  were  randomly  selected.  Table  rv  shows 
a comparison  of  the  results  of  these  trajectories  for  the 
rigid  surface  and  for  the  moving  surface  cases.  It  is 
not  clear  why  the  moving  surface  atoms  seem  to  increase 
the  total  probability  of  molecule  formation  by  a factor  of 
2.  The  energy  transferred  between  the  gas  atoms  and 
the  surface  atoms  for  both  the  equal  and  unequal  mass 
cases  was,  in  the  mean,  4%-7%  of  the  total  energy  avail- 
able to  the  product  molecules.  The  mean  distribution  of 
the  available  energy  among  vibration,  rotation  and  trans- 
lation Is  changed  very  little  by  the  motion  of  the  surface. 
More  work  is  needed  to  assess  carefully  the  role  of  sur- 
face motion  on  angular  distributions  (which  measures 
momentum  transfer  and  are  frequently  a more  sensitive 
probe  of  collision  dynamics)  as  well  as  the  interaction 
of  surface  atoms  and  adsorbate  atoms  with  comparable 
masses.  The  generalized  Langevln  oscillator  model  of 
Adelman  and  Doll”  could  be  Included  to  mimic  the  ef- 
fects of  the  full  lattice  on  the  reaction  zone.  The  pres- 
ent studies  indicate  that  the  product  state  distributions 
of  recombining  atoms  are  not  very  sensitive  to  the  mo- 
tion of  the  surface  atoms. 

VI.  CONCLUSION 

From  the  above  calculations  we  conclude  that  many  of 
the  rules  postulated  by  Folanyl  for  gas-phase  energy  re- 
quirements and  disposal  mechanisms  are  applicable  also 
to  recombination  of  adsorbed  atoms  for  gas-solid  inter- 
actions. We  have  found  that  repulsive  potential  surfaces 
which  in  the  gas-phase  require  vibrational  energy  In  the 
bond  under  attack  for  reaction  do  not  lead  to  reaction  for 
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heterogeneous  recombination.  This  is  true  also  for  en- 
dothermic potential  surfaces.  Attractive  potential  sur- 
faces give  rise  to  more  vibrational  excitation  in  the 
product  molecules  and  less  translational  energy  than  do 
mixed  energy  release  surfaces.  However,  the  light- 
heavy  mass  combination  was  found  to  remove  the  in- 
crease in  vibrational  excitation  found  in  attractive  rather 
than  mixed  energy  release  surfaces.  The  energy  is 
channeled  instead  into  rotational  excitation.  The  mbced 
energy'  release  surfaces  give  rise  to  more  translational 
energy,  consistent  with  Polanyi's  result.  The  rigid  sur- 
face was  found  to  be  a good  approximation  for  calculating 
product  state  distributions  in  these  recombination  studies. 
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Introduction 


There  have  been  experimental  indications  that  for  atomic  recom- 

1 2 

bination  on  solid  surfaces,  energy  is  disposed  of  in  highly  specific  ways-  ’ 

For  example,  in  recombination  of  oxygen  on  a variety  of  metal  surfaces,  the 

data  indicated  that  most  of  the  reaction  energy  was  often  deposited  into  the 

2 

gas  phase  diatomic  molecule.  This  was  in  contrast  to  conventional  ideas  that 
there  was  essentially  equilibrium  partitioning  of  energy  between  gas  and  solid 
for  atomic  recombination,  Recent  data  by  Hal  pern  and  Rosner^  clearly  confirm 
the  specific  nature  of  energy  disposal  in  heterogeneous  reactions  for  recombi- 
nation of  atomic  nitrogen  on  a variety  of  metals.  The  data  is  sufficiently 
precise  and  reproducible  that  some  hypothesis  about  the  chemical  dynamics  of 
recombination  reactions  on  surfaces  could  be  proposed.  It  had  been  known  for 
some  time  that  when  molecules  react  in  the  gas  phase,  various  modes  of 
energy  disposal  (i.e.,  as  vibration,  rotation,  or  translation  in  the  products) 
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are  not  all  equivalent.  Whether  or  not  a large  fraction  of  the  reaction  energy 
appears  as  product  translation  or  vibration  is  known  to  be  determined  by  the 
detailed  energy  contour  along  the  reaction  path,  relative  masses  of  the  atoms, 

4 

and  other  subtle  features  of  the  potential.  A reaction  occurring  on  a solid 
surface  shares  many  of  the  same  features.  A reaction  path  can  be  constructed, 
the  energetics  along  this  path  can  be  described  and  the  dynamics  of  the 
resulting  reaction  can,  in  principal,  be  studied.  The  sole  complicating 
feature  is  the  presence  of  a solid  surface  upon  which  the  reaction  occurs. 
Therefore,  it  is  not  surprising  that  heterogeneous  reactions  would  exhibit  pre- 
ferred modes  of  energy  disposal.  Precisely  how  the  solid  surface  complicates 


2 


the  dynamics  of  heterogeneous  reactions  is  far  from  clear,  but  the  fact  that 
interesting  dynamical  effects  exist  in  heterogeneous  reactions  seems  beyond 
dispute. 

One  of  the  most  useful  consequences  of  specific  energy  disposal  in 

gas  phase  reactions  has  been  the  formation  of  vibrational 1y  excited  product 

molecules  and  resulting  infra-red  chemical  lasers.^  The  question  naturally 

arises  whether  reactions  promoted  by  a heterogeneous  catalyst  can  lead  to 

equally  useful  vibrational  population  inversions.  Various  experimental 
1 2 

groups  * have  been  led  to  speculate  that  a large  fraction  of  the  available 
reaction  energy  does  appear  as  product  vibration,  but  the  experimental  support 
is  extremely  fragile. 

In  analogy  with  similar  studies  for  gas-phase  dynamics,  we  have 
undertaken  a series  of  computational  studies  attempting  to  model  vibrational 
excitation  in  heterogeneous  reactions.  Our  method  for  approaching  this 
problem  has  followed  closely  the  theoretical  methods  (i.e.,  classical  trajec- 
tories on  a series  of  model  potentials)  whereby  the  behavior  of  gas  phase 
chemical  lasers  was  elucidated.  Sinc^  these  methods  were  productive  for  gas 
phase  chemical  lasers,  there  is  reason  to  believe  they  will  also  lead  to 
important  insights  into  heterogeneous  chemical  dynamics  and  vibrational  exci- 
tation. Because  of  this  close  analogy,  we  think  it  is  productive  to  review 
briefly  the  history  of  gas  phase  chemical  lasers  to  point  out  the  reasons  for 
our  current  studies,  and  by  learning  from  history,  we  hope  to  indicate  likely 
future  directions  for  heterogeneous  chemical  dynamics. 

Polanyi  and  co-workers  observed  vibrational  population  inversion 
via  infra-red  chemiluminescence  in  HCl  formed  from  the  reaction  of  K + Cl2- 
The  possibility  of  constructing  an  infra-red  laser  was  explicitly  discussed 
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in  a paper  submitted  for  publication  in  I960.®  An  extensive  series  of  classi- 
cal trajectory  studies  on  a model  potential  surface  (LEPS)  was  undertaken  in 
order  to  understand  the  conditions  under  which  vibrational  population  in- 
version could  be  expected  in  exothermic  reactions.  The  first  results  were 
submitted  for  publication  in  1962,^  and  a detailed  discussion  of  energy  re- 
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lease  in  exothermic  reactions  appeared  in  1966.  Various  features  in  the 
potential  energy  surface  in  addition  to  exothermicity  were  shown  to  be 
necessary  for  efficient  population  inversions  to  occur.  Acting  on  Polanyi's 

O 

1960  prediction,  Kasper  and  Pimentel  produced  the  first  chemical  laser  in 
1965  using  the  H2/CI2  reaction.  The  function  of  a heterogeneous  catalyst 
is  to  alter  the  energy  profile  along  the  reaction  path  to  accelerate  the  rate 
of  formation  of  products.  For  gas  phase  reactants  and  products  the  overall 
exothermicity  cannot  be  affected  by  the  catalyst.  But  such  catalysts  alter 
reaction  rates  drastically  by  altering  the  energetics  of  intermediate  reaction 
states.  Therefore,  from  atom-diatom  collisions  in  the  gas  phase  we  know  the 
importance  of  details  in  the  energy  path  for  determining  final  state  popu- 
lations and  can  ask,  how  will  a catalyst  affect  these  final  state  population 
distributions?  Clearly,  there  must  be  an  effect  since  the  energy  profile 
along  the  reaction  contour  changes,  but  what  can  be  expected  in  detail? 
Following  what  was  dona  for  gas  phase  chemical  lasers,  we  chose  to  use  classi- 
cal trajectories  to  attempt  to  elucidate  these  features  in  the  gas-solid  inter- 
action, most  likely  to  lead  to  population  inversions.  Don  Bunker  was  one  of 
the  pioneers  in  classical  trajectory  methods  and  would  not  be  surprised  to 
find  the  method  applicable,  essentially  unchanged,  to  the  emerging  area  of 
dynamics  of  heterogeneous  catalysis. 


4 


In  this  paper,  we  present  some  recent  results  concerning  chemical 
dynamics  at  the  gas-solid  interface  via  classical  trajectories.  This  repre- 
sents the  culmination  of  a series  of  papers  whose  goal  was  to  assess  the 
possibility  of  obtaining  vibrational  population  inversion  by  means  of  a 
heterogeneous  chemical  reaction.  Therefore,  in  addition  to  reporting  our 
recent  results  concerning  the  dynamics  of  adsorption,  we  will  provide  an 
overview  and  rationale  for  the  recent  work  in  this  area. 

In  order  to  model  realistic  chemical  dynamics,  a reasonable  model 
potential  is  needed.  We  used  a modified  4-atom  London-Eyring-Polanyi-Sato 

(LEPS)  model  potential  to  account  for  the  interaction  of  a diatomic  molecule 

9 10 

with  a solid  surface.  In  our  first  studies,  the  surface  was  assumed  to  be 

rigid  and  provide  simply  a substrate  upon  which  the  reaction  could  occur.  This 

is  consistent  with  the  simplest  picture  of  a catalyst  as  providing  merely  a 

lower-energy  pathway  from  reactants  to  products  but  (as  a rigid  surface)  not 

adding  or  removing  energy  from  the  reaction.  A very  similar  model  potential 

was  developed  Independently  by  Gelb  and  Cardillo^^  and  used  to  investigate 

12 

isotope  effects  for  hydrogen,  and  cluster  formations  on  a copper  surface. 

Consistent  with  the  rigid-surface  model,  we  have  assumed  that  the 
reaction  energy  to  be  partitioned  in  the  product  molecules  must  come  entirely 
from  the  reactants  rather  than  from  the  substrate.  For  efficient  population 
inversion  to  occur,  considerable  energy  must  be  available.  Since  catalysts 
are  typically  heavy  metals,  one  does  not  expect  a very  efficient  heat  transfer 
between  the  surface  and  rapidly  reacting  molecules  on  the  surface  as  an  inter- 
mediate state  in  a fast  reaction.  Therefore,  the  source  of  the  reaction  exo- 
thermicity  is  assumed  to  be  the  reaction  itself  and  the  surface  plays  a minor 
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role  dynamically  (but  a major  role  in  altering  the  energetics  of  the  reaction 

path).  This  assumption  has  been  tested  numerically  and  found  to  be  accurate 
13 

to  'lO*.  As  discussed  above,  experimental  results  have  appeared  in  which 

the  fractional  energy  deposited  in  the  solid  has  been  measured.  While  the 

experiments  are  not  precisely  comparable  with  the  present  calculations,  the 

data  clearly  shows  that  strong  chemical  bonds  are  frequently  formed  without 

1 2 

significant  energy  transfer  to  the  solid.  ’ 

Section  II  provides  a description  of  our  model  for  chemical  dynamics 
at  the  gas-solid  interface,  emphasizing  the  possibility  of  obtaining  vibra- 
tional population  inversions  from  such  a system.  Section  III  discusses  our 
most  recent  results  for  the  dynamics  of  adsorption  as  related  to  the  possibili- 
ties of  obtaining  vibrational  population  inversions  and  connects  these  with 
previous  results  for  the  dynamics  of  heterogeneous  recombination.  Conclusions 
are  drawn  and  the  possibility  of  heterogeneously  catalyzed  reactions  leading 
to  useful  population  inversions  are  discussed  in  Section  IV. 
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II. Description  of  the  Model  , 


We  envision  studying  a reaction  of  the  general  form 


I 

i 


A2  + 02 


(gas  phase) 

» 2AB(V,J) 


> 2AB(V’,J') 

(catalyst) 


(la) 

(lb) 


In  practice,  the  direct  reaction  (la)  frequently  has  a substan- 
tial activation  barrier  and  correspondingly  slow  rate.  The  catalyst  is 
chosen  specifically  to  remedy  this  situation  but,  by  altering  the  energy 
contours  between  reactants  and  products  the  (V,J)  distributions  will  also 
change.  Conceptually,  it  is  convenient  to  discuss  the  catalytic  reaction  (1b) 
ir  two  steps.  In  the  first  step,  a diatomic  molecule  (A2  or  B2)  collides 
with  the  solid  surface  and  dissociatively  chemisorbs.  In  the  second  step,  the 
adsorbed  atoms  A and  B encounter  one  another  on  the  surface  and  desorb  as  AB. 

We  do  not  consider  the  case  in  which  a gas  phase  species  adsorbs  as  an  un- 
dissociated diatomic  molecule  and  this  species,  rather  than  chemisorbed  atoms, 
participate  in  the  reaction.  Nor  do  we  consider  direct  reactions  between  gas 
phase  A2  and  adsorbed  B atoms  (or  the  inverse).  Both  of  these  cases  (par- 
ticularly the  second)  could  be  important  for  practical  catalytic  systems  oper- 
ating at  high  pressure.  However,  the  additional  complications  introduced  into 
the  reaction  path  by  having  to  follow  more  atoms  materially  increase  the  compu- 
tational effort.  Since  our  understanding  of  the  dynamics  is  rather  crude,  we 
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believe  It  is  justified  to  study  the  simplest  realistic  system.  Hence,  we 
consider  reaction  (1b)  to  proceed  via  dissociative  chemisorption  of  both  A 
and  B,  a subsequent  encounter  of  the  adatoms  on  the  surface,  and  finally, 
the  desorption  of  AB.  Such  a model  would  be  appropriate  at  low  pressures 
and  low  coverages  of  adsorbed  species  and  (obviously)  for  a system  in  which 
A2  and  Bg  do  dissociatively  chemisorb.  In  summary,  our  2-step  process  is 

A2.B2  (gas)  ^ A,B  (chemisorbed  atoms)  >-AB(V,J),  (gas)  (2) 

Step  I Step  II 

The  requirement  that  dissociative  chemisorption  takes  place  means 
Step  I should  be  exothermic  for  both  A2  and  B2.  As  argued  above,  we  require 
AB  to  receive  its  internal  energy  from  the  overall  reaction  exothermicity, 
rather  than  from  the  catalyst.  Therefore,  we  have  the  further  restriction 
that  the  overall  reaction  I + II  be  exothermic.  The  detailed  procedure  for 
modeling  the  overall  reactions  I and  II  is  as  follows: 

Step  I:  A single  diatomic  molecule  (either  A2  or  B2,  for  simplicity 
we  will  assume  A2)  collides  with  the  surface  of  the  solid.  Energy  may  be  ex- 
changed with  the  solid,  the  A2  bond  may  break,  or  the  molecule  may  rebound 
into  the  gas  phase.  This  calculation  Involves  Just  one  diatomic  molecule 
colliding  with  a moving  surface  and  is  well  within  the  computer  technology 
developed  previously.  This  part  of  the  calculation  provides  estimates  of: 
the  efficiency  with  which  dissociative  chemisorption  occurs,  the  distribution 
of  kinetic  energies  for  the  atoms  adsorbed  on  the  surface,  and  the  relative 
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fraction  of  energy  transferred  into  the  solid  (and  thus  lost  to  the  reaction). 
All  are  calculated  for  a variety  of  potential  surfaces  constructed  consistent 
with  our  assumptions  (i.e.,  exothermic  for  St*p  I). 

Step  II:  Two  adsorbed  atoms  A.B  encounter  each  other  on  the  sur- 
face, recombine  and  from  a gas-phase  diatomic  molecule.  Again,  we  must  con- 
sider only  two  atoms  interacting  with  a moving  surface  to  compute  the  state 
distribution  of  the  products  and  the  energy  lost  to  the  solid.  The  only 
direct  connection  between  Step  I and  Step  II  is  through  the  distribution  of 
kinetic  energies  of  adsorbed  atoms  on  the  surface.  However,  judging  from  gas 
phase  dynamics,  the  formation  of  vibrational ly  excited  product  molecules  seems 
to  be  more  sensitive  to  the  energy  profile  of  the  reaction  contour  in  the 
formation  step  rather  than  to  the  detailed  initial  conditions  of  the  reactants. 
That  is,  for  an  exothermic  reaction  A + (B,C)  the  product  state  distribution 
is  largely  determined  by  the  energy  profile  along  the  reaction  path,  and 

effects  due  to  the  initial  state  of  the  (B,C)  or  due  to  the  kinetic  energy  of 

14 

collision  are  not  as  important. 

To  this  approximation.  Steps  I and  II  can  be  decoupled  and  computed 
separately.  Detailed 
ready  appeared^ and 


studies  of  recombination  dynamics  (Step  II)  have  al- 
will  be  briefly  summarized  below.  In  this  paper  we 
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it 
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report  computational  results  for  the  adsorption  step  (I). 


III. 


Results  and  Discussion;  Dynamics 


As  indicated,  there  were  reasons  to  feel  that  Step  II,  the  recom- 
bination step,  would  be  the  most  crucial  in  determining  the  population  dis- 
tribution of  the  final  molecules.  Therefore,  this  step  was  studied  first  and 
in  some  detail. We  briefly  summarize  the  results  here.  Using  a model 
potential  based  on  the  London-Eyring-Polanyi-Sato  (LEPS)  method,  nine  potential 
surfaces  were  generated  (Figures  1,  2).  The  original  rigid-surface  model  was 
generalized  to  allow  for  motion  of  the  surface  atoms,  and,  hence,  gas-solid 
energy  exchange.  Three  of  the  9 surfaces  were  endothermic  for  recombination 
(Number  I,  II,  III),  while  the  remaining  6 were  exothermic  for  recombination. 
Various  barriers  and  modes  of  energy  release  (early,  late)  were  included. 

The  results  for  the  9 surfaces  are  summarized  in  Tables  I and  II.  The  3 endo- 
thermic surfaces  did  not  lead  to  efficient  recombination  nor  to  significant 
reaction  energy  appearing  in  vibration  even  though  kinetic  energy  considerably 
in  excess  of  the  minimum  was  supplied.  Hence,  one  requirement  for  efficient 
energy  deposition  in  AB  seems  to  be:  although  A2,  82  dissociatively  chemisorb 
on  the  material,  AB  should  not  dissociatively  chemisorb  (i.e.,  AB  recombination 
in  the  surface  should  be  exothermic).  This  should  be  a useful  guide  in 
choosing  candidate  systems  for  experimental  studies.  We  see  that  some  of  the 
systems  studied  (i.e.,  VII),  produce  vibrational  population  inversions,  quite 
comparable  to  the  HF  laser.  Also,  system  VII  has  a satisfactory  amount  of 
energy  in  V*2,  which  is  another  important  characteristic.  Unfortunately,  the 
best  results  were  obtained  for  the  case  of  equal  mass  which  is  likely  to  be  a 
homonuclear  molecule  not  suitable  for  IR  emissions  (although  CO  could  approxi- 
mately qualify). 
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In  summary,  the  results  for  the  recombination  seem  to  indicate  a 
good  population  inversion  may  be  obtained  if:  (1)  the  recombination  is  exo- 
thermic, (2)  atoms  have  roughly  equal  masses,  and  (3)  from  potential  VII, 
the  energy  is  released  early  in  the  recombination  with  no  barrier  to  recombi- 
nation. That  vibrational  excitation  is  produced  more  readily  by  an  attractive 
potential  (i.e.,  VII),  is  consistent  with  the  work  of  Polanyi^^  on  gas  phase 
A + BC  collisions. 

It  remains  to  examine  the  adsorption  step  to  see  if  significant 
energy  loss  to  the  surface  occurs.  We  have  shown  in  previous  studies  that  ad- 
sorption will  proceed  efficiently  even  without  energy  transferred  to  the  solid 

10c 

if  it  is  exothermic  and  does  not  encounter  an  activation  barrier.  Also, 
we  do  not  expect  the  detailed  motion  of  the  atoms  A and  B across  the  surface 
to  be  as  important  as  the  potential  in  determining  the  final  state  distribution 
of  AB  product  molecules.  Therefore,  our  chief  concern  is  with  the  energy 
transferred  to  the  solid  during  adsorption  and  whether  or  not  such  energy 
transfer  could  be  a serious  loss  of  energy  otherwise  available  for  AB  exci- 
tation. We  have  previously  studied  the  energy  transfer  for  H2,  and  HD  colliding 
with  a tungsten  surface  using  a potential  that  had  a barrier  to  adsorption. 

Here,  we  report  results  for  a few  potentials  lacking  such  activation  barriers, 
as  well  as  for  some  heavier  particles  with  early  and  late  energy  release, 
and  for  heavier  gas  molecules.  Potentials  1-3  in  Figures  1 and  2 are  exo- 
thermic to  adsorption  (endothermic  to  recombination),  provide  both  early  and 
late  energy  release,  and  these  were  used  for  our  study.  Both  light  and  heavy 
mass  combinations  were  used  (to  simulate  H2  and  a system  with  the  mass  of  Br2 
colliding  with  tungsten)  and  the  results  are  reported  to  Table  III.  Both 
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Isolated  diatomic  molecules  were  described  by  a Morse  curve  with  potential 

a 

parameters  given  previously.  This  potential  was  chosen  to  reproduce  the 
properties  of  the  isolated  H2  molecule  and  only  the  mass  was  increased  to  80 
amu  to  mimic  a heavy  molecule  colliding  with  the  surface.  Hence,  the  system 
labeled  "Br2"  has  only  the  mass  of  "Br"  and  the  other  properties  of  dissoci- 
ation energy,  force  constant,  and  equilibrium  separation  appropriate  for  H2. 

The  Morse  parameters  for  the  isolated  diatomic  molecule  affect  the  LEPS 
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potential  in  a complicated,  non-linear  way.  Hence,  attempting  to  change 
the  Morse  parameters  to  model  more  closely  a physical  Br2  molecule,  would 
change  the  total  potential  and  the  energetics  of  the  reaction.  We  believe 
it  is  more  important  to  produce  a comparison  on  the  same  potential  with 
simply  a heavier  molecule. 

Both  vibrational  and  rotation  were  Included.  Approximately  1.75  eV 
of  vibrational  energy  is  present  in  both  systems  (H2(V=3),  "Br2"(V=31 )). 

Similarly,  the  rotational  quantum  numbers  were  selected  such  that  equivalent 
rotational  energy  (0.09  eV)  was  present  for  both  molecules.  For  all  cases 
studied,  9 surface  atoms  were  permitted  to  oscillate  about  their  lattice  sites, 
according  to  the  Einstein  model  discussed  previously. The  parameters  for 
this  oscillation  about  the  lattice  are  the  same  as  those  used  in  the  recombi- 
nation study  of  reference  13a. 

Ejqj  in  Table  III  is  the  total  energy  available  to  be  depositied  into 
the  solid;  adsorption  energies  plus  kinetic  and  internal  energies  of  the  gas 
molecule..  The  surface  temperature  was  initially  300°K.  The  energy  gained  by 
the  solid  because  of  this  adsorption  is  AE^.  It  is  calculated  by  summing  the 
kinetic  and  potential  energies  of  the  solid  atoms  as  they  move  about  their  lattice 
sites.  This  is  not  a precisely  defined  number  because  the  couplings  between  the 
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the  surface  atoms  and  adsorbed  atoms  will  cause  energy  to  flow  back  and 
forth  until  the  adsorbed  atoms  recombine  and  leave  the  surface,  or  the  ad- 
sorbed atoms  leave  the  region  of  moving  surface  atoms.  In  this  case,  the 
trajectories  were  followed  until  the  final  separation  of  the  adsorbed  atoms 
exceeded  the  region  of  moving  surface  atoms. 

From  Table  III  we  see  that  the  larger  mass  will  increase  the  energy 
transfer  by  roughly  a factor  of  2-3  although  the  mass  increases  by  a factor 
of  80.  For  the  worst  cases,  the  energy  transfer  is  approximately  10%  of  the 
total  energy.  However,  there  are  recent  indications  that  the  Einstein  model 
used  ht?re  tends  to  underestimate  energy  transferred  to  the  solid  when  compared 
with  generalized  Langevin  (GL)  computations.^^  The  GL  computations  considered 
only  one  moving  surface  atom,  and  it  is  not  clear  how  much  this  affected  the 
results.  That  is,  several  Einstein  oscillators  may  provide  a closer  approxi- 
mation to  a more  accurate  GL  model  (there  being  more  degrees  of  freedom  in  the 
solid  in  which  to  store  energy)  than  a single  oscillator  would  provide  to  a single 
oscillator  GL  model.  Until  these  calculations  are  performed,  the  results  must 
remain  somewhat  uncertain.  To  be  conservative,  AE^/E^qj  in  Table  III  could 
be  too  small  by  a factor  of  2.  For  ^2  + W(OOl)  the  net  energy  loss  to  the 
solid  is  still  within  the  original  estimate  of  10%.  However,  for  the  mass  80 
system,  20-25%  energy  loss  could  occur. 
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Conclusion 


In  conclusion,  the  current  series  of  computer  simulations  indicates 
that  vibrational  excitation  should  be  possible  for  heterogeneous  reactions. 

The  conditions  seem  to  be  that  both  steps  in  the  reaction,  dissociative  chemi- 
sorption and  recombination,  should  be  exothermic.  Efficient  population  in- 
version occurs  for  early  energy  release  in  the  recombination  step  and  for 
nearly  equal  mass  atoms.  A question  that  has  not  been  fully  explored  is  the 
energy  loss  to  the  solid  for  all  conditions  of  interest.  The  present  compu- 
tations indicate  that  this  could  be  a serious  loss  (i.e.,~20%)  for  heavy 
systems,  but  should  be  quite  tolerable  for  light  gases.  There  are  indications 
that  heating  the  substrate  would  tend  to  reduce  the  net  energy  loss  into  the 
solidJ® 

One  is  tempted  to  speculate  about  the  possibility  of  constructing 
a ":atalytic  chemical  laser"  given  the  estimates  of  population  inversions  given 
above.  However,  this  endeavor  is  not  likely  to  be  fruitful  at  this  stage  be- 
cause there  if  simply  too  little  data,  either  experimental  or  theoretical, 
with  which  to  do  a realistic  estimate.  The  present  research  indicates  that  a 
vibrational  population  inversion  should  be  possible  in  reactions  catalyzed 
by  a solid  surface.  There  are  also  indications^®  that  collisions  of  vibrationally 

excited  molecules  with  surfaces  do  not  invariably  lead  to  deactivation.  De- 

-4 

activation  probabilities  as  low  as  10  are  not  uncommon.  This  would  be  the 
chief  loss  mechanism  not  already  present  in  chemical  lasers,  and  it  may  be 
possible  to  overcome  it  by  an  appropriate  design.  The  overriding  need  now  in 
heterogeneous  chemical  dynamics  appears  to  be  for  good  measurements  of  final 
state  population  distributions  for  a variety  of  prototype  systems.  If  classical 
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trajectory  studies  such  as  these  in  any  way  stimulate  interest  in  such 
measurements,  or  provide  guidance  as  to  how  they  should  be  performed,  the 
work  and  career  of  Don  Bunker  in  developing  the  technology  must  receive  due 
credit. 
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FIGURE  CAPTIONS 


Figure  1 - Equipotential  contours  (in  eV)  for  the  approach  of  H2  to  W(001) 
for  the  different  potential  surfaces  studied.  The  contour  levels 
are  equally  spaced  for  any  given  potential  surface,  X indicates 
the  position  of  the  reactants,  two  atoms  adsorbed  in  adjacent  ICN 
sites  (5CN  for  IX).  • locates  our  choice  of  the  "transition  state." 
For  details  of  the  construction  of  these  surfaces,  see  Ref.  13a. 

Figure  2 - Energy  profiles  along  the  reaction  paths  for  the  potential  surfaces 
given  in  Figure  1.  The  dashed  line  marks  the  region  of  transition 
from  reagents  to  products  (denoted  by  • in  Figure  1).  (From 
Reference  13a.) 
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TABLE  I.  (From  Reference  13a.)  Mean 

distribution  of  energy  in  product 
molecules:  (a)  rigid  surface; 

(b)  two  surface  atoms  free  to  move. 


Potential 

% Molecule 

* ^ib 

* ^ot  * 

c 

% Energy 

formation 

trans 

Transfer 

(a)  Equal  masses 

IV 

(a) 

41 

41 

15 

44 

(b) 

72 

49 

11 

34 

6. 32 

V 

(a) 

36 

60 

20 

20 

(b) 

78 

52 

26 

15 

6.69 

VI 

(a) 

42 

51 

10 

39 

(b) 

65 

57 

8 

30 

5.24 

VIII 

(a) 

30 

68 

18 

14 

(b) 

69 

62 

18 

13 

7.04 

(b)  Light-heavy  mass 

combination 

IV 

(a) 

21 

47 

13 

40 

(b) 

73 

51 

10 

34 

4.79 

V 

(a) 

40 

49 

30 

21 

(b) 

79 

37 

40 

19 

3.  86 

VI 

(a) 

40 

51 

13 

36 

(b) 

69 

47 

17 

32 

5.89 

VIII 

(a) 

26 

37 

40 

23 

(b) 

70 

46 

27 

23 

4.26 
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TABLE  II.  (From  Reference  13a.)  Product 
vibrational  state  distribution 
{%)  summed  over  rotational 
sublevels. 


Potential  v * 0 1 2 3 ‘ 


(a)  Equal  masses 


IV 

18 

39 

30 

12 

1 

0 

V 

11 

17 

36 

34 

3 

0 

VI 

8 

36 

39 

15 

2 

0 

VIII 

1 

20 

48 

26 

5 

0 

IX 

14 

29 

32 

24 

1 

0 

(b)  Light-heavy  mass  combination 


IV 

11 

20 

20 

30 

15 

4 

V 

15 

21 

26 

22 

7 

9 

VI 

io 

10 

23 

23 

25 

10 

VIII 

18 

23 

13 

13 

27 

8 

IX 

39 

22 

14 

20 

3 

2 

21 


TABLE  III.  Energy  lost  to  the  solid  for 
potentials  I,  II,  III  defined 
in  Figures  1,  2. 


Potential  Surface 

'rel'*''! 

V 

J 

^^s^^TOT  {%) 

Hg  + W(OOl)  I 

0.3 

0 

0 

3.2 

I 

0.3 

3 

0 

1.6 

I 

CO 

o 

0 

3 

2.9 

II 

0.3 

0 

0 

4.5 

II 

0.3 

3 

0 

2.3 

II 

0.3 

0 

3 

4.3 

III 

0.9 

0 

0 

3.0 

III 

0.9 

3 

0 

2.3 

III 

0.9 

0 

3 

2.9 

III 

0.3 

3 

0 

1.9 

III 

0.5 

3 

0 

2.4 

"Br2"+W(001)  I 

0.3 

0 

0 

6.2 

I 

0.3 

31 

0 

5.0 

I 

0.3 

0 

30 

7.0 

II 

0.3 

0 

0 

10.2 

II 

0.3 

31 

0 

13.1 

II 

0.3 

0 

30 

10.6 

III 

0.9 

0 

0 

6.5 

III 

0.9 

31 

0 

10.3 

III 

0.9 

0 

30 

8.4 
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ABSTRACT 


The  generalized  Langevin  method  for  gas-solid  scattering  is 
combined  with  a LONDON-EYRING-POLAMYI-SATO  type  potential  function 
in  order  to  allow  calculation  of  collisions  and  simple  reactions 
at  the  gas-solid  interface.  Ener^  transfer  and  dissociative  ad- 
sorption for  selected  collision  configurations  for  a diatomic  mole- 
cule scattering  off  a Cu  (100)  single  crystal  face  have  been  studied 
by  integrating  the  classical  trajectories.  Special  emphasis  is 
placed  on  determining  the  relative  importance  of  the  microscopic 
parameters  of  the  collision  system  on  the  microscopic  mechanisms  of 
these  processes.  Energy  accommodation  and  dissociative  adsorption 
were  found  to  be  greatly  effected  by  the  Debye  temperature  of  the 
solid. 
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I.  IflTRODUCTIOri 

Previously,  the  many-bo4y  nature  of  gas /sol Id  processes  has  prevented  ] 

the  theorist  from  developing  a realistic  dynamical  picture  of  gas/solld  I 

collisions.^  The  recent  development  of  the  generalized  Langevin  method  for  i 

gas/solld  processes  reduces  this  many-body  problem  to  a coirutatlonally 

2 3 4 

practical  few  body  problem.  ' * The  generalized  Langevin  theory  Is  based 

( 

on  a few  key  Ideas.  Since  the  Impinging  gas  molecule  Interacts  strongly 
with  only  a few  of  the  solid  atoms,  the  dynamics  of  these  solid  atoms  must 
be  explicitly  accounted  for.  For  the  collision  processes  of  Interest,  the  > 

motion  of  the  remainder  of  the  lattice  Is  observable  only  by  Its  Influence 

on  the  region  of  strong  Interaction.  Therefore,  one  need  not  explicitly 
follow  the  detailed  motion  of  the  surrounding  lattice  as  long  as  a suffi- 
ciently large  region  of  strong  Interaction  has  been  Included.  The  surrounding 
lattice  ente^  the  dynamics  via  friction  kernels  and  random  forces.  The 
many-body  problem  also  manifests  Itself  In  the  diatom-solid  surface  potential 
energy  problem.  This  problem  has  been  cast  Into  a computationally  tractable  > 

form  by  the  use  of  a London-Erying-Polanyl-Sato  type  potential  function  developed 
by  McCreery  and  Wolken.®  Combining  this  potential  function  method  with  the 
generalized  Langevin  method,  we  present  calculations  of  collisions  and  sinple 
reactions  at  the  gas/solld  Interface  In  which  both  chemically  realistic  model 
potentials  and  a treatment  of  the  full  lattice  dynamics  has  been  Included.  | 

As  a simple  example  of  a chemical  reaction,  we  choose  the  dissociative  ad- 
sorption of  a diatomic  molecule  on  the  (100)  single  crystal  face  of  copper. 

Also,  the  energy  transfer  occurring  during  the  scattering  of  diatomic  mole- 
cule for  a solid  Is  Investigated  via  collision  configurations  arbitrarily 
constrained  to  be  non-reactIve.  ? 
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This  paper  is  Intended  to  be  a preTirainary  study  of  the  dynamics  of  di- 
atomic mo1ecule/so11d  reactive  processes  with  particular  attention  to  the 
lattice  effects.  Therefore,  we  have  approximated  the  directly  struck  zone 
(primary  zone)  of  the  surface  by  one  moving  solid  atom.  This  solid  atom  was 
further  confined  to  oscillate  perpendicularly  to  the  surface  plane.  Also, 
the  center  of  mass  of  the  diatom  perpendicularly  approaches  the  moving  solid 
atom  In  all  our  trajectories.  Only  broadside  and  perpendicular  collision 
configurations  were  used  for  the  approach  of  the  bond  of  the  diatom.  Such 
configurations  greatly  simplify  the  elucidation  of  the  microscopic  mechanisms 
Involved. 

In  this  study,  we  wish  to  gain  a better  understanding  of  how  the  micro-  ' 
scopic  parameters  of  the  collision  system  affect  the  microscopic  mechanisms 
of  gas/sbild  processes  in  which  chemical  reactions  occur.  In  the  scattering 
of  atoms  from  solid  surfaces,  the  duration  of  the  collision  Is  known  to  deter- 
mine the  extent  of  participation  of  many  body  dynamics  in  energy  acconmodation.^ 
At  a given  beam  energy,  the  more  massive  the  gas  atoms  the  smaller  the  velocity 
and  the  longer  the  collision  time.  In  the  present  study,  we  are  able  to  vary 
collision  times  by  using  fictitiously  large  gas  atom  masses.  The  mass  ratio 
of  gas  atom  to  surface  atom  also  affects  the  extent  of  energy  exchange.  Para- 
meters of  the  solid  such  as  the  surface  temperature  and  the  Debye  temperature 
of  the  solid  have  known  effects  In  gas  atom/surface  energy  accommodation. 

The  relative  Importance  of  the  parameters  of  the  collision  system  on  diatom/ 
surface  collision  process  Is  poorly  understood  due  to  the  possibilities  of 
dissociative  chemisorption  and  the  participation  of  Internal  molecular  modes. 

The  present  study  addresses  these  problems. 
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In  Section  II,  we  describe  the  generalized  Langevin  equation  method  along 
with  the  specifics  of  the  potential  energy  surface.  In  Section  III,  we  discuss 
our  calculations  and  results. 


II.  GLE  METHOD 


The  generalized  Langevin  theory  for  gas /sol  Id  processes  and  Its  numerical 

2 3 4 

Implementation  Is  extensively  discussed  elsewhere.  * * Here,  we  briefly 
review  only  those  aspects  of  the  theory  pertinent  to  the  specifics  of  our 
calculations.  The  set  of  generalized  Langevin  equations  (GLE)  Is 

^j(t)  +1  B^j(t  - T)rj(T)dT  + 


"J'j 


R'j(t)  “ WCR'^  »R2*r'^  ,r2». . . ) 


(1) 


Tj  denotes  collectively  the  coordinates  of  the  atom  j In  the  solid  (j  ■ 1, 
2,...)  while  R^(1  ■ 1,2}  denotes  the  coordinates  of  the  two  atoms  comprising 
the  Incident  gas  molecule.  The  present  calculation  Is  for  a diatomic  mole- 
cule Incident  on  a solid  surface,  but  this  Is  not  an  Inherent  restriction  In 
the  theory.  Is  the  danplng  kernel  accounting  for  the  dissipation  of 
energy  from  solid  atom  J to -the  remainder  of  the  solid.  Rj(t)  Is  the  random 
force  experienced  by  solid  atom  j due  to  Interactions  with  the  surrounding 
solid.  In  general,  m^  Rj  could  be  different  for  different  atoms  In 
the  solid,  but  this  generality  Is  not  used  In  the  present  calculations  as 
we  are  attempting  to  model  the  collisions  of  a diatomic  gas  with  a solid  of 
uniform  composition.  Therefore,  the  subscripts  j will  be  omitted.  For  the 
present  calculation,  the  damping  kernel  e(t)  Is  approximated  by  a one-term 
damped  sine  function  obtained  by  a numerical  fit  to  a modified  bulk  Debye 
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model  of  the  solid  which  was  shown  previously  to  be  a reasonable  approxima- 


tion/*^ 


Uq  Is  the  fundamental  harmonic  frequency  of  the  solid,  also  within 


solid 


(2) 


the  modified  Debye  approximation.  The  random  force  R(t)  Is  assumed  to  be 

Q 

gausslan  and  generated  by  procedures  described  elsewhere. 

W Is  the  potential  between  the  Incident  diatomic  molecule  and  the  solid 
surface.  It  was  chosen  to  have  the  form 

W - W4+  I (W 
movin 
solid 
atoms 

The  Interaction  of  the  diatomic  molecule  with  a rigid  solid  surface  (frozen 

lattice)  Is  taken  to  be  a modified  4 Body  LONDON-EYRING-POLANYI-SATO  (LEPS) 

5a 

type  potential  function  W^.  To  allow  for  motion  of  the  solid  atoms,  we 

add  W|y|_^,  a 3-Boc|y  LEPS  function  between  the  moving  solid  and  the  diatom. 

To  avoid  overcounting  the  Interaction  of  the  molecule  with  the  rigid  surface, 

we  then  subtract  the  value  of  the  3-Body  LEPS  at  the  equilibrium  position  of 

the  solid  atom,  Wu  « , and  sum  over  all  moving  surface  atoms  analogous  to 

”"^eq 

a procedure  described  In  detail  elsewhere  where  was  modified  by  pairwise 
additive  potentials  In  contrast  to  the  3-Body  correction  used  here.®^ 

The  application  of  LEPS  potential  functions  to  the  Interaction  of  mole- 
cular hydrogen  with  the  (100)  face  of  Cu  has  been  described  previously  by  >. 


Gelb  and  Cardlllo  In  the  form: 
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The  hydrogen  molecule  interactions  and  are  taken  to  be  tte  usual 
Mor^e  and  anti-Morse  functions  respectively,  and  the  hydrogen  atom-surface 
interactions  are  given  by 
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with 
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and 


■ ^cunf  - f - -rj  “«(^) 

The  hydrogen-copper  potential  parameters:  Sj^^,  ^-Cu*  Cu* 

Z^uH»  Y^t  and  ; are  given  in  tables  I and  II  of  reference  9.  This 

potential  has  a minimum  energy  for  dissociation  of  7.0  Kcal/mole  and  a barrier 
to  diffusion  of  10  Kcal/mole.  The  hydrogen- copper  potential  parameters  for 
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the  3>Body  LEPS  potential  were  chosen  to  be  those  of  Gregory  et  al^®,  and  they 
are  given  In  table  I.  . 


III.  CALCUUTIONS  AND  RESULTS 

Energy  exchange  and  dissociative  adsorption  have  been  studied  for  a di- 
atomic molecule  colliding  with  a solid  surface  by  Integrating  classical  tra- 
jectories for  the  potential  function  described  here*  approximating  H2  Cu(lOO). 
Since  our  primary  Interest  Is  In  understanding  what  microscopic  parameters 
most  sensitively  effect  gas/solld  collision  processes,  we  vary  these  parameters 
to  elucidate  their  effect. 

(A)  ENERGY  EXCHANGE 

Energy* exchange  with  the  surface  was  studied  by  colliding  a diatonic 
molecule  with  the  solid  restricting  the  bond  axis  of  the  molecule  to  be  per- 
pendicular to  the  plane  of  the  surface,  and  positioning  the  Incoming  molecule 
directly  over  the  struck  solid  atom.  For  diatom  Intemuclear  distances  close 
to  the  equillbruua  distance,  this  collision  configuration  Is  repulsive,  and 
thus  no  dissociative  adsorption  should  occur  for  the  energies  studied  here. 

The  purpose  of  restricting  the  collision  configuration  Is  to  remove  several 
complications  such  as  multiple  collisions  with  the  surface  and  diatom  rota- 
tional dynamics.  The  Internal  energy  of  the  diatom  was  chosen  to  be  the  same 
as  that  of  the  ground  vibrational  and  rotational  state  of  the  hydrogen  mole- 
cule, and  the  kinetic  energy  of  the  diatoms  center  of  mass  was  selected  to 
be  10  Kcal/mole.^^  This  value  was  selected  to  simulate  the  experimental 
conditions  of  ialooch  et  al.®  The  temperature  of  the  solid  In  their  experi- 
ments was  850"K.  Sets  of  100  trajectories  were  run  for  each  set  of  solid/ 
beam  conditions. 


1 

I 


i 


The  energy  loss  of  a beam  of  structureless  gas  particles  colliding  with 

4 

a solid  is  expected  to  decrease  with  increasing  solid  tenperature.  This  was 
conducted  for  a molecule  with  internal  vibrational  structure  by  running  sets 
of  100  trajectories  for  the  solid  at  temperatures  of  0*K,  300*K,  and-  850* K. 

The  diatom  was  a fictitious  hydrogen  molecule  which  was  composed  of  atoms 
20  times  as  massive  as  a hydrogen  atom.  We  shall  call  this  molecule  H(20). 

The  results  are  shown  in  Fig.'  1,  and  as  expected,  average  energy  transfer 
does  decrease  with  increasing  solid  temperature.  It  is  very  interesting  to 
note  that  these  collisions  were  vibrationally  adiabatic;  that  is,  there  was 
no  loss  of  vibrational  energy  by  the  diatomic  molecule.  This  seems  to  be 
the  result  of  the  large  amount  of  zero  point  vibrational  energy  present  in 
the  hydrogen.  The  extent  of  the  many  bocjy  effect  on  energy  exchange  may  be 
determined  by  removing  the  damping  kernel  and  random  force  from  the  equations 
of  motion.  This  model  is  called  the  single  oscillator  or  "Einstein"  model 
and  the  results  are  also  shown  in  Fig.  1.  The  average  energy  exchange  is 
decreased  by  a factor  of  roughly  2 although  the  trend  is  remarkably  well  dup- 
licated. 

If  one  varies  the  Debye  temperature  of  the  solid,  one  effectively  varies 
the  primary  spring  constant  and  the  scaled  time  for  energy  dissipation  into 
the  solid.  This  is  due  to  the  scaling  of  the  constants  in  the  damping  kernel 
by  the  Debye  frequency  and  that  the  primary  frequency  ^ wq.  One  ob- 

4 

tains  Uq  from  the  Debye  temperature  of  the  solid.  As  one  Increases  the  Debye 
temperature  one  increases  the  stiffness  of  the  spring  which  should  decrease 
the  average  amount  of  energy  transferred  from  the  molecular  beam.  However, 
an  increased  Debye  temperature  also  means  a quicker  dissipation  of  energy 
into  the  solid.  We  have  tested  the  effect  of  increasing  the  Debye  temperature 
of  the  solid  by  running  sets  of  trajectories  for  the  model  H2-CU  system  where 
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the  surface  Debye  temperature  was  set  at  100**K,  236'’K»  and  350"K.  Again, 
the  diatom  was  H(20).  The  solid's  temperature  was  300'’K.  The  results  of 
our  calculations,  see  Fig.  2,  shov/  that  the  main  effect  of  Increasing  the 
Debye  temperature  Is  to  decrease  the  average  energy  transfer.  The  single 
oscillator  model  results  for  the  average  energy  transfer  were  again  too 
small  but  with  the  correct  trends. 

The  general  trend  of  Increasing  gas  mass  Is  expected  to  Increase  the 
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average  energy  transfer  for  cube  models.  The  opposite  trend  would  be 

13 

predicted  by  Goodman's  driven  oscillator  models.  For  a 300*K  solid  with 
copper's  Debye  temperature,  we  studied  the  energy  transfer  of  the  hydrogen 
molecule,  and  fictitious  molecules  with  atoms  10  times  and  20  times  the 
mass  of  hydrogen  in  H(10)  and  H(20).  He  observe  that  the  average  energy 
transfer  shov/s  a peak  for  the  fictitious  molecule  H(10),  see  Fig.  3.  As 
expected,  the  hydrogen  molecule  sho^/s  the  least  amount  of  energy  transfer 
In  the  6LE  results.  The  decline  In  energy  transfer  for  more  massive  gas 
molecules  Is  the  result  of  the  struck  solid  atom  recoiling  and,  upon  rebound, 
depositing  energy  again  Into  the  slower  moving,  more  massive  molecules.  The 
single  oscillator  model  should  converge  to  the  GLE  results  as  one  reduces 
the  collision  time,  I.e.  reduces  the  mass  of  the  diatom  at  constant  kinetic 
energy.  In  Fig.  3,  we  see  that  the  GLE  and  single  oscillator  results  for 

i 

the  hydrogen  molecule  are  very  similar.  Only  the  effect  of  the  random  force 
causes  the  slight  difference  In  average  energy  transfer. 

(B)  DISSOCIATIVE  ADSORPTION 

Ue  also  studied  the  dissociative  adsorption  of  a diatomic  molecule  on 
tlie  (100)  single  crystal  face  of  copper  as  a function  of  solid  temperature, 
Debye  temperature  of  the  solid,  and  gas  mass.  Again,  sets  of  100  trajectories 
v/ere  run  for  each  of  the  different  solld/beam  conditicr:.  To  simplify  the 
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slmulatlon  of  adsorption,  we  restrict  the  configuration  of  the  diatom  to  be 
parallel  to  the  surface  with  the  center  of  mass  perpendicularly  approaching 
the  oscillating  solid  atom.  The  angle  a line  parallel  to  the  surface  makes 
with  the  bond  axis  of  the  diatom  was  Monte. Carlo  averaged  with  the  other 
Initial  conditions.  This  configuration  permits  the  diatom  to  dissociate  upon. 

Impact  with  the  surface  with  Its  constituent  atoms  becoming  trapped  In  wells 
of  atomic  adsorption.  The  Importance  of  lattice  dynamics  on  dissociative 
adsorption  becomes  Immediately  obvious  In  this  configuration.  This  Is  easily 
seen  by  the  fact  that  the  rigid  surface  appt*ox1mat1on  results  lii  no  absorp- 
tion for  any  of  the  three  cases  of  gas  mass  studied  with  their  Internal  energy  ^ 

equivalent  to  the  ground  vibrational  and  rotational  state  of  hydrogen. 

Experimentally  It  has  been  observed  that  changes  In  the  temperature  of  i 

the  solid  will  have  the  largest  effect  on  dissociative  adsorption  when  the  < 

dissociation  Is  proceeded  by  physisorptlon. Our  potential  surface  does  not 
have  a well  for  .physl sorption  and  therefore  we  expect  little  temperature  effect 
on  our  results.  The  effect  of  temperature  and  dissociative  adsorption  was 
studied  using  the  fictitious  molecule  H(20).  The  GLE  curve  shows  little  var- 
iation with  temperature,  see  Fig.  4.  However,  the  single  oscillator  model 
results  are  marked  by  the  small  probability  of  adsorption  at  0‘K.  A sample 
study  of  the  trajectories  Indicates- that  this  Is  due  to  spurious  processes. 

The  diatom  actually  dissociates  and  the  atoms  bounce  off  adjacent  solid  atoms 

and  recombine.  Since  the  single  oscillator  model  contains  ^l^of  ttie-energy 

lost  by  the  colliding  molecule,  ^e  solid  atom  Is  able  to  kick  out  the  recom- 
bined gas  molecule.  This  mechanism  occurs  In  the  GLE  calculations  when  the 
gas  atoms  bounce  back  before  the  surface  solid  atom  can  dump  Its  energy  Into 
the  lattice.  Either  quickly  moving  gas  atoms  such  as  hydrogen  or  long  Deb/e 
periods  of  the  solid  bring  this  about.  Therefore,  one  sees  the  Importance 
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many  bo4y  effects  have  on  sinple  0*K  dissociative  adsorption.  Another  spurious 
mechanism  occurs  for  finite  temperature,  single  oscillator  trajectories.  The 
atoms  of  the  dissociated  gas  molecule  oscillate  on  the  surface  In  potential 
wells  adjacent  to  the  solid  atom  with  which  they  collided.  The  moving  solid 
atom  kicks  some  of  these  oscillating  gas  atoms  Into  adjacent  wells.  Our  model 
contains  no  energy  loss  mechanism  to  the  solid  when  the  adsorbed  atoms  are  at 
such  a great  distance  from  the  moving  solid  atom.  None  of  the  6LE  trajectories 
tested  showed  this  spurious  effect. 

Increasing  the  Debye  tenperature  of  the  solid  has  two  counterbalancing 
effects.  As  one  Increases  the  Debye  tenperature  one  Increases  the  spring 
stiffness  and  therefore  decreases  the  average  amount  of  compression  of  the  metal 
atom  oscillator  for  the  same  energy  Input.  When  the  extent  of  compression  of 
the  solid  atom  Is  decreased,  the  attraction  of  the  atoms  of  the  colliding  diatom 
for  atomic  adsorption  wells  is  decreased.  Also,  the  atoms  of  the  colliding 
molecule  will  have  more  kinetic  energy  on  the  average  after  collision  with  a 
stiffer  spring,  as  was  shown  by  the  ener^  transfer  calculations  (Fig.  2).  The 
rate  of  energy  dissipation  Into  the  solid  Increases  with  Increasing  Debye  tem- 
perature. 

In  order  to  ascertain  the  combined  effect  of  these  factors,  we  ran  reactive 
trajectories  for  several  Debyt  tenperatures  for  the  fictitious  molecule  H(20). 

Our  calculations  Indicate  that  Increasing  the  Deby^  te^erature -of  the  -solid 

— decreasis'  the"  whbimt  of  dissociative  adsorption,  see  Fig.  5.  It  Is  Interest- 
ing to  note  that  several  of  the  GLE  trajectories  for  the  Debye  tenperatures  of 
100*K  and  350*  were  able  to  recombine  and  desorb.  The  single  oscillator  model 
displays  no  dissociative  adsorption  for  the  Debye  temperature  of  350*IC.  With- 
out the  many  body  dynamics,  the  single  oscillator  model  quickly  becomes  equi- 
valent to  the  rigid  surface  resulv^. 
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As  Fig.  6 displays,  an  Increase  In  the  gas  mass  results  in  a steady  in- 
crease in  adsorption  probability  for  both  the  6LE  and  single  oscillator  model. 


SUriMARY  AND  CONCLUSION 

We  have  presented  calculations  of  collisions  and  simple  reactions  at  the 
gas/solid  interface.  As  a preliminary  study  of  the  dynamics,  we  restricted 
the  collision  configurations  to  remove  certain  dynamical  complications.  In 
this  light,  we  have  studied  the  relative  importance  of  the  microscopic  para- 
meters of  the  collision  system  on  the  microscopic  mechanisms  of  dissociation 
and  energy  accommodation.  A large  dependence  of  both  energy  accommodation  and 
dissociative  adsorption  on  the  Debye  temperature  of  the  solid  was  found.  An 
Increase  In  the  Debye  temperature  drastically  decreases  both  the  energy  trans- 
fer and  the  amount  of  dissociative  adsorption.  Energy  accommodation  was  found 
to  decrease  with  increasing  solid  temperature;  however,  dissociative  adsorp- 
tion, which  Is  not  proceeded  by  physl sorption.  Is  not  affected  by  this  Increase. 
The  Importance  of  Including  the  many  body  nature  of  the  solid  was  shown  through 
the  spurious  results  of  the  single  oscillator  model  In  the  dissociative  ad- 
sorption studies. 

In  conclusion,  we  wish  td^MphasIzif  that  this  study  Is  otrly  a start  In  — 
unraveling  the  microscopic  mechanisms  of  reactive  gas/solid  processes.  In  a 
more  realistic  collision  model,  several  surface  atoms  will  be  strongly  Inter- 
acting with  the  Incident  gas.  A further  complication  of  such  a study  would 
be  the  determination  of  the  number  of  moving  surface  atoms  which  must  be 
included.  We  have  started  our  studies  with  a more  simplified  model  so  as 
not  to  obscure  the  nature  of  the  process. 
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TABLE 

The  parameters  of  the  Morse  functions  which  are  used  in  the  3 
BODY  LEPS  for  the  gas  atonncopper  atom  potential  are  taken  from 
A.  R.  Gregory,  A.  Gelb,  and  R.  SilbeyJ®  We  modified  their  more 
parameter  associated  with  the  well  width  in  order  to  prevent  the 
tunneling  of  gas  atoms  into  the  solid. 

Well  Depth  • 0.316  eV 
Well  Width  -0.5  a.u. 

Equilibrium  Intemuclear. Distance  • 2.3  A 
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FIGURE  CAPTIONS 

% 

AE  vs.  temperature  of  the  solid  for  H(20)/Cu(100)  scattering. 

-1/2 

Error  estimates  are  the  usual  error  in  the  mean  N ' (standard 
derivation)  where  N is  the  number  of  trajectories. 

AE  vs.  0JJ  for  H(20)/Cu(1000)  scattering. 

AE  vs.  gas  mass  for  di atom/ Cu(  100)  scattering.  The  diatoms  are 
hydrogen,  H(10),  and  H(20). 

Sticking  probability  vs.  the  temperature  of  the  solid  for 
H(20)/Cu(100)  scattering. 

Sticking  probability  vs.  0p  for  H(20)/Cu(100)  scattering. 

Sticking  probability  vs.  gas  mass.  The  diatoms  are  hydrogen, 
H(10),  and  H(20). 
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^ Figure  2. 

^ Fi gure  3. 

1 

I Figure  4. 

r 

j Figure  5. 

1 Figure  6. 


so 


so 


